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2 PROFESSOR K. PEARSON ON THE INFLUENCE OF NATURAL 


(1.) Introductory. On the Influence of Selection upon Correlation. 


AT an earlier stage in the development of the statistical theory of evolution it 
was suggested that the coefficient of correlation (GAuron’s function) might be found 
constant for all races of the same species—in fact, it was considered possible that this 
coefficient might be the long-sought-for criterion of identity in species. Professor 
We pon, following up this suggestion of Mr. Ganron’s, then made the elaborate 
series of measurements on crabs with which his name will always be closely associated. 
To a first approximation these researches seemed to confirm the possibility of 
GALTONS function being a true criterion of species. When, however, a finer 
mathematical test was applied to Professor WeLbowN's observations as well as to 
other statistical series for organs in man,” it became clear that the coefficient of 
correlation varied from local race to local race, and could not be used as a criterion 
of species. A slight investigation undertaken in the summer of 1896 convinced me 
that the coefficient of correlation between any two organs, is just as much peculiar and 
characteristic of a local race as the means and variations of those organs. In fact, 
if local races be the outcome of natural selection, then their coefficients of correlation 
must in general differ. The object of the present paper is to show, not only that 
natural selection must determine the amount of correlation, but that 1t 1s probably 
the chief factor in the production of correlation. If selection, natural or artificial, 
be capable of producing correlation, then 1t seems impossible to regard all correlation 
as evidence of a causal nexus, although the converse proposition that all causal 
nexus denotes correlation, is undoubtedly the most philosophical method of regarding 
causality. 

Tn dealing with the influence of selection on correlation, I shall suppose the distri- 
bution of complex groups of organs to follow the normal correlation surface—the 
generalised Gaussian law of frequency. I shall further assume the selection surfaces 
to be normal in character. Neither of these assumptions is absolutely true, but the 
Gaussian law in a good many cases describes the frequency sufficiently closely to 
enable us to obtain fair numerical results by its application. Probably in all cases it 
will enable us to reach qualitative if not accurate quantitative theoretical deductions. 
I have the less hesitation in asserting this, as Mr, G. U. Yure has recently succeeded 
in deducing the chief formule for correlation and regression as given by the 
Gaussian law from general principles, which make no appeal to a special law of 
frequency.{ 

* “Mathematical Contributions to the Theory of Evolution.—11.,” ‘Phil. Trans., A, vol. 189, pp. 266 
and 280. See also E. WARREN, “Variation in Portunus depurator? ‘Roy. Soc. Proc.” vol. 60, 
pp. 233-4. 

+ See a series of letters in ‘ Nature,’ vol. 54, 1896, arising from a discussion upon a paper by A. R. 
WALLACE, 

t “Roy. Soc, Proc., vol. 60, p. £77. 
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In particular, if a selected group be not given with very great accuracy by a 
normal frequency surface, still we may, I think, consider ourselves justified in 
supposing that the effects of the actual selection, and those of a normal selection 
with the same means, the same amounts of variation, and with correlations of the 
same intensity will be at least qualitatively alike in character, if they be not indeed 
exactly the same quantitatively within the limits fixed by the probable errors of the 
constants. 

My plan in this memoir will be as follows :—I shall first state the fundamental 
theorem in multiple correlation with a new proof, so that the formulee required may 





be once for all collected for reference.* I shall then give the algebraic investigation 
of the new formule for selection. I shall afterwards consider special simple cases, 
and illustrate them by examples. Finally I shall draw attention to the nature of 
the selective death-rate as indicated in cases of this kind, and consider at length its 
algebraic theory. Throughout I shall endeavour to illustrate the somewhat complex 
algebra by arithmetical examples. 


(2.) On the Fundamental Theorem in Multiple Correlation. 


I have shown in my memoir on “ Regression, Heredity, and Panmixia” (‘ Phil. 
Trans. A, vol. 187, p. 261) that if the n variables of a complex be functions of m 
(m > n) independent variables with frequency distributions following the normal law, 
and such that the principle of superposition holds for the deviations from the means 
supposed small; then the frequeney of the complex with deviations from the 


means of the n variables lying between x, and a, + 8x,, %, and a, + Sa,... a, 
and £, + dx, will be z ôx, dv)... Sx, where z: 

SE 1 9 y i 

z = zy ekpt. — ¿lO (opt) + 2B) o o . a . DÄ 


Here 2, Ga Cpg. - - are constants, and S, denotes a summation for every value of p, 
and X, for every pair of values of p and q in the series from 1 to n. 

In the same memoir (p. 302) I have determined the values of Zos Cop Cpg U terms of 
the correlations 7,, and the standard deviations o, and o, of the n variables. This 
pomt had already been considered by Professor Eperworra (‘ Phil. Mag., vol. 34, 
p. 201, 1892), and some further results by Mr. A. BLACK, reached before his death in 
1893, were published in the ‘Camb. Phil. Trans.’ (vol. 16, p. 219, 1897). The 
present investigation is, I think, novel, and adds to results already reached others 
required in the present memoir, so that it thus places together with a fairly simple 
proof all the fundamental results to which I shall have occasion to appeal later. 


* We have used these formule for several years, but they do not appear to have been hitherto published 
in a collected form. 
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Let us consider the quadric of the n™ order 
(==, (Ct) + 28, (ee) == constant, 


and fix our attention on two of the variables, say the first two, zy and x, If these 
be considered constants, the quadric of the remaining n — 2 variables will not now 


be referred to its “centre.” But its centre, tz, £}... %, will be given by the 
equations : 


Cry, + Cagla F Cass lr. AE Ca = O. 


CUL F Cala F Cala F Gout, Cy = O. 


Cinti +F Cono SS Cans T Canta T : : : zb Canen = 0. è S S (11.). 
The following expressions will not be zero, but will be written « and £ : 

CUE] F Cala F Cala F Cala ho... Oy = a. 

Cioli T Cola F Cys TS A il Bo... (111.). 
Now, if A be the discriminant, 


Ci, Ci, C135 ‘ : : Cin 
Cor, Coo, Cos, e e e Con 


C31, C395 C33, ° R . Coy : V a a è : d (iv. ) H 


d 





Got 3 az: Cis e ` ° Cx 2 


and C,, the minor corresponding to c,,, we have by solving the above n linear 


equations, 


sa 











A EE 
: A 2 A 
whence: Ca eS eS A A CY. 
Dun — Oig” p 1 Oy — Le 
Generally also : 
Eo Gree SI T PCa 
Ly — A 





> Kl y E 
(Cy Loa — Cag@q9)@ + (Cy, Cag — CiO) ta 
Cy Cog — Cro” 


— D x + Dy Ae, Say, e e e D H e D D D $ e (v1.). 
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This determines the central co-ordinate for any variable x, for a given value of 


x, and Zo. 


Now let us transfer the quadric to £a, £4.. 


Q = Sal Ly Le: SR Coogl -{- Gaata +. 

= Sia" (e, H xy) (ue F Cogo F Calg H 

EE Lo g f 
i S + Catat e 
+ 2 (Oo + Cala E Cat +... 

] 

\ + ya's... 
+ a, ES + Calo F Clg +... 
i F Cpl tH 


Making use of the n linear equations (11.) and (11) : 


Q = 5723 (£y + 07) (Cals + Cta t- 
+ x (a + egs + Cyt, +... 
+ 2o (B + Cog F Cka H 


ge / / / 
= BH Uy (Cag 3 Ar yyy Tt 


+ ax, + Bar, 


. L, 48 OFISTN, 


It may be written 


, Lo, 


+ (aa A 


T Geh 


y 


SR Crta ) 
n 
\ 


2 / 
"E: Ci 10 


+ Cra 


F e 


F GE 
+ CNA a) 


+ nt. 


) 
+ Cnn) 


For arranging vertical columns in rows, the remaining terms are 


al: Loan + Cao + SLES (Cogg) 
Lal, fou + Cata + DI Lg 


Si Le (Et, A Cuak F SEN (Cngq)), 


each line of which vanishes by the equations (11.) for the centre, 


Accordingly : DECH ‘a 
= Xy Los 


where Q’ is a quadratic function of £g, Cy . . 
Hence 2 = Bee ee 


(vil. ). 


. x’, not involving x, and x, at all. 


Now integrate z with respect to all the variables o w,,... a", from — œ to -+ œ, 


keeping xv, and x, constant. 
Then, although the origin is a function of x, and x,, 


+O fo +0 ee i ; A 
| | DEER CORO a UN 
— un 


— 0 — 00 


cannot involve x, and æ, but only gau Cag.» «Eng» &C. ; let the result be Zo 


E EE Es Bra 


Then : 
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must be the correlation surface for «, and %, for all values of £s, mu... Uy. We may 
write it in the form 


2 = e expt. =r 4 A BEE EE S LS l : (viii). 
( 11 C22 Gd Cr 





Comparing this with the known form of the correlation surface for two variables, * 








3 t 1 1 dy i 2119119 dd E ; 
SE = a ER A Een ee REENEN (ix.) 
210,0 1 a CEP SS ml Led 9 ; 


ERT ës oy 
we have at once 


a es C 11 Cog — Un fal ae Cy Coy = Uy 
Sone PRA G 


2 Fa 
On (1 Tig AC ACh > 


O A A O A 
Whence ieee ln A a ee lee 


=D A, Giese Ane, ee oe an ee Eh 
Or, generally : oe 8 
m OVE Cr) in ep oe et RR ee EECH 


SNE 0 ETE A 


Thus correlations and variations are fully determined in terms of the discriminant 


abd its minor for the comstants en Cag, Cross Cops Cap Cpgs -> Cow 


We have next the inverse proposition to find the co in terms of the 1% and o's. 
. We have, by well-known propositions in the theory of determinants : 


4 y y 
Got + CeCe e Cirat... BE Ge, 
N 
Gun F CC + CCa . . E CnC = 0. 
i Y { 
CaCa + CrO + Gag, E 2 . +. Oz = 0. 


nta + CEs + EE Sp oes zb ST Ee N 


o 8 
Or, Gu F Cpe? 12010 + Cpu zk, zb: Ceuta = E 
5 
211,100) Se Ci zk Cig oT Ae oo a SS EE on == 0; 


S 4 2 z al 
Cfa A CiaP 090% F Cogs zk, . . A Cin. andan = 0 


CM F Cy yay, zk C.na zb . + A EE 
Hence solving 


CLIO, o A 


* PEARSON, ‘ Phil. 'Crans.,' A, vol. 187, p. 264. 
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| 2 Ki E] 4 ; 
where SE On, TiTa, 130103, ee Padda ` 
e | 
T210 102, Oz, Taatfotfo, 6 6 e Von Oy, 
o | 
M3190 1O3, Tautfotfz Ts, — 6 Tania, 
| | 
| 
a na ch 2 
| 1 vi Ojon, ) 7120 CT y, 1 213738 > . ` . T; 





and 5, 18 the minor corresponding to the constituent r,,07,0,. 


But EE EE e 


el oes SE 2 
yy ses OF To Gr sn a Ty sn 6 Oy, Rugs 


where R is the determinant, 





| | H Vy, WÉI Pin | 
Tan z 123) Pon 
| May Ia L; Pon | 3 
| a H (xvi. i 
a D ) 
H Ai y H OI GER s S K ii 





and R,, is the minor corresponding to the constituent r,,. Thus we have 


1 N E S 
C 11 < 9 le Cig — ee. $ R/R . . ° . . (vi, ), 
Si EE 


or, generally, 
Fo EE SE RaR (xvii) 
4 PG GRO, PY D D a D D a o Po 


bole 


: E ? N 1 HI L =p pra "4 y Dé A D 
2 = Zu expt. — S, (=> 2) + 28, ae E (x1x.). 
R R Ou R Tn Ty 


b 
\ 7 
N D 


It remains to determine zu from the fact that the volume of the surface = N. 


GL o +n Lo ; 
| 
Or, N = % | | SES | expt. (— $Q) du, ds... dx, 
dm y mm — Y 
O ppo +0 j "H E No) 
Se EE e hh 
mee Zoff fa e . e O; | | e e e | expt. GE 2 di D Lp -+ 255 R LE yt r) d 
Ye dD — 0 \ Au 


E ZS SE 
EE e Et, Eh 
by writing top = Ga Se. 
Now, integrate first with regard to al, writing first 


* + Phil, Trans. A, vol. 197, p: 302. 
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(Lo, Hi 
s (Pe e +28, (0) 


“0 / Cal pe Pan J a 
al o (ke 29 + 5 ghee ct R Cn 
11 
> 2.3 fT. a T : y 
+ > Se II ap! 2) + 23 Rogby — acht, gl op! 


where 2, is the sum for all values of p, and 3, for all pairs of values of p and q, from 
2 to n inclusive, 


H = e n a (Lap 
= e, + A E + 28 (0 aya!) where E, = 3, ex sl 


Zu 





For if R = Ry, or the determinant of the correlation coefficients, omitting all 
involving the first variable, 2.e., the first row and column of R, the dca eee Ro 
corresponding to the minor of tho constituent Zon m HI or to the second minor (. Ri)» 


of R is given by 


Se RB’, = (RyRy, — BARBAR, (xxi), 
and tor p = q, | 
A Gan) 


Hence, integrating x’, between the limits -+ œ and — œ, we have 


CH 
N — Zoff soe MS or e 
Lo Pen +0 i 7 Gen 5 Rng rar’ j 
4 ) ¢ oy f oe SE eat 
x | | e o o y expt. SE ES È ( R >/ SI +- KH JI R 20) AL ALE o e e de DÉI 
= OO — OD — OO f 


This is of precisely the same form as before, except that we have the factor 
/ Ze vi R/R’, and the multiple integral is reduced by one integration and by the 
disappearance of all correlations involving the first variable. Now, integrate with 


regard to al, The sole effect will be to multiply by a factor V 2r /R//R”, where 
R” is the minor of R’ not involving correlations of the second variable. Thus, by 





repeating the process, we have ultimately 





N=>200,...0, (Vr VRIR VRIR” / BTR”... ¿Je n 


N 
(= a == oy 
"7 erg eeng al R 


or, (xxiii. ), 


which gives the constant of the surface. 
The preceding investigation enables us also to deal with two further points. 
(a) Given n variables, what is the mean value m, + x, of the first variable and 


its variability o”, for definite values m, + ho, Ma + ha... Ma + h, of the other 
(n — 1) variables ? 
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Clearly, when «’,«’,... x’, are constants, the distribution of a: is of the form 


J Ry, Së A 
constant X expt. — y R a. + aN = Lp ; 


Or, re-introducing the ex, ga... c, we have a distribution about the point given by 


H ry? ` 
ar Ryo kg R, 3 Jig ie Ly \ 
A ee as. O a ERA as a 3 = E be SE 
Zb e 1 D r | a a e R J H 
a EE 1 On 


(xxtv.). 


with standard deviation 


(b) Given n variables, what are the mean values m, + 2, My, + Zo, the standard 


deviations o”,, o”,, and the correlation 1,, of two of them, when we give definite 
l> 2 19 , 
values Mz + hz, Mi + hi... . Mm, +h, to the remaining (n — 2) variables ? 
In this case we have from (1.) 


= expt. — $ [0,2 + 2c 42% Xa + Coro? 
+2 (cils + Culi EL... + Goal 
+f 2 (tags + Cay Hb... Lehm: 


+ terms not involving æ; anda}... . (xxv.). 


Writing K, for the coetticient of x,, and K, for that of x,, we have for the centre 

















a (Kye, = Kaca) OR Wes 08 bo + Ciy K3) 
e EE Q ? EE C 2 g 
Deag "Eug ou 22 Dua 
ep 
de a > e d \ 
ES Circa ~ Dote SC Lo = — Y E h 
SE Ee SE SS ere Fp? 
1122 2 11622 12 
3 A + P 4 
E 5 (Ra — RoRo hy ARA > Roplty — GE hy 
a? E An i 
, Hg — Mig” Fy ty Fog — Ey o 


by transferring to the minors of R and the os. Or, 





gen, Zen , 
2, = — 012, KE = C= = EE E er DEE NT 
Here R” is the determinant formed by striking out the first two rows and 
columns of R; p”,, is the minor obtained by striking out the second row and column 
from R, and then the first row and mi column; py, the minor obtained by striking 
out the first row and column, and then the second row and p column. But a 
comparison with (xxiv.) shows us that these values for a, and i, are precisely what 
we should have obtained for the regression equations of the lst and 2nd variables 
respectively alone on the other n — 2 variables. ~Thus the existence and the 
correlations of x, have no effect on the value of x,, nor those of x, on the value of æ, * 
VOL. CC.—A. e c | 
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Returning to (xxv.), we remark that the terms of the second order in x, and x,, on 
which the correlation and variations depend, are not altered by a transfer to the 
centre x, and x, of the array. 

Hence by (ix.) and (xvu.) we have 


= Uu The. EE +s 
ie == <A os A ee e (RANA. 
13 degt E ( ) 


This is the partial correlation of the 1st and 2nd organs for the remaining n — 2 
organs with constant values. 
Again, 
1 1 


Ca SAA lo = III] áÁA 
E E Eg (ae) 
Whence we easily find from Een and (xxvi) 


e 3 D 
R Rag cap ao 3 R Ra 


KN ff 2 2 O viel ` os O es a a 5 
DH E i 
j ` RaRa — ui 


CV SO: > Doa? 
Ru, ES Rye 





Or, 


Si, = 0, Ball, 7. = oy vi HAN . . . . (xxviii), 


where R” is, as before, the determinant R without its first two rows and columns. 


These by (xxiv.), are what we should have reached by ignoring «e, in finding ol, and 
x in finding o”, 
(3.) General Theorem in Selection. 


To find the selected means, the selected varvitions and selected correlations, when 
y organs are selected, naturally or artificially, out of a complex of n organs. 

Let the selected group of q organs have their means raised h,, hy, hs, .. . hy (some 
of these quantities may be negative); their standard deviations changed from 
o a... 07 tO Sp ën 83,... Sy and their mutual correlations from ms 743, .. + Tip 
Togs Togs ++ Tas, + + Pg LO Pro Piss ++ + Pio Pas Paw +++ Paps + Parar 

The whole system of n organs before selection will be defined by the means as 
origin of measurement for each organ, by the standard deviations oe, Uy, Og, .. . Oy, 
and by the coefficients of correlation 7, fis, Me EE A IR E 
R be the determinant 
| en ioe, Vises ee ee, Oar 


| Pa, a do a de Ree 


t 


a nr Wi 
| Pais Tex Pus > E, l 





and R,, the minor corresponding to the constituent r,e Then the unselected 


population is given by the frequency surface of equation (x1x.), 
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Ray dë Ro Cabo ; 
z = constant X expt. — $ 18; E ES Gi + 28’, E a vet | pee E 


A 
D R opr; 


For brevity, we can also write this in the form 


z = constant X expt. — 4 {S} (60,7) + 28) (092) . `, (xxix. is 


Now consider for the time only q + 1 organs—namely, the first q organs and the 
organ (m > q), and let us write R (u), if u be > q, for the determinant : 
Rus e E EENE E SE, AS): 
Tols 1 > Pox, : i : Por Va, 


> Ki 3 gr 
Par H 325 1 > : : R 1 Zu? Pu 


Vou Tyr, Tase à 1, Ton 
Tuis Pou Puss S : Pug 1 


Then if R (u)pp be the minor corresponding to the constituent Zus, and if 
p — Bike 1 
CN E aaa pr 


frequency distribution 


the distribution of the q +1 organs will be given by the 


2 = constant X expt. — 4 {S (by) + 2S, (Opp Bp Gyr). (XXX). 


S, being a sum for every value of p’ throughout the q + 1 organs, and S, for every 
pair of values. 
Now let the first g organs be given values h,, ha, ... Ma then the mean value of æ, 
will be given by 
/ 


A A Don 
SEENEN 


UU uu 


R (u) u SE R u 2u Cu R U Jon u A 
om - Le l ht: Lol P. hy t.. SEET, 


Ra o Bisi 0, Be (xxxii). 


Now these coefficients can be found at once if q be known. 
For example : 














SSC E R (De Ge 
3 R (Wer lau > 
g=, — Bh — De tate RO — Za "Tir 
, = == = a 
E (un l-re’ SE (U Juu Le 
q=3 = R (u) e __ w AU — Peg) — Piou E E Pos (Pieu F TigPo) 
; = EE 
Jh L — Ti — Pag — Pig + 2779799713 


R R (2 )ow Ta R Dik 


and 
E CH SS R (2 EK 


can be written down by symmetry .  (xxxlil). 
C 2 
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Now suppose in the expression (xxix.)™ we were to put £u, Za, ... «e, equal to 
hi, hy, . . . hy respectively, then the system of equations to find the means of 
Carro Lagos + «+ Corn, © + + La for this array would be the n -- q equations 

a. dÉ S T DE 
Cr gri la F Co, 941 fto F ++ Hey gi he F Corr, Ca F gaara P ooo Kage Cn = O, 


l Gë / ¿A 
C1, 949 a F Co 949 to F > > F Cg gag Mg Cr, nun F Co, quo Cgo ho > + EC quo Un = 9. 


9 


Cin d SS Con hs SE SE F Cy, i h, T Cg+1,n are it Co 42,0 euro E + Cn xl = 0 ° (XXXIV. ), 


where Cour Cura... 2, are the “co-ordinates of the centre” of the array of the 
n — g organs. 

If we were to solve these equations we ought to get precisely the solution for 
e, (u>q<n+ 1) that we have found for a”, m (xxx1.) above, where none of the 
coefficients involve correlation-coefficients other than those of the first q organs 
among themselves and with the adi organ. ‘This result flows from the pretty obvious 
law that the mean of the ah organ for an array determined by values of the 
first g organs cannot be in any way dependent on our considering the relation of this 
selection of q organs to any additional organs beside the u'”: see p. 9. 

Thus the solution of (xxxiv.) is simply obtained by putting u = q + 1, q + 2, 
q + 3,... successively in (xxx11.). 

Let us now select the first q organs not with absolute values, but varying about 
means fy, Aa... A, with standard deviations e, $3... Ssp and with mutual corre- 
JALIONS pn pia tte Para e e pra) Veevhavecthien tommultyply Eeer 
by an exponential quadratic function of the 1, +/h,, v7 -+-hy,... 2,+h,, t.e., the selective 
correlation surface, and divide it by another exponential quadratic surface, 1.e., the 
primary correlation surface of the q organs %,, #,...«, This follows, since the 
frequency of each complex of n organs must be reduced in the ratio of the selected 
to the primary frequency of the complex of q selected organs. But it will be clear 
that such a reduction must give us a result of the following form for the final 


frequency surface of the n organs : 


Z = constant X expt. — kloe + Cot? +... + Cto 
SE GE EEN SE eg 3° SR Gef 
+ Bëss Lk... (Vv and gd <q+1) 
ck Dep, ... (0<q+1and u > 9) 
+ Zeyl bk... (w and u” > 9) 


+ linear terms in £y, Ze (v and <q +1) . (xxxv.), 


where the c’s denote the changed cs, 
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Now if we differentiate the quadric to find its “ centre” we have n equations in 
Lio. Ly... Ly, but the solutions of these, 1f £+; . . £a were eliminated, are known 
to be the “centre” hi, hy... Agp Hence we require only n — q equations involving 
dann... L, and we can put the As for the remaining values æ... æ, Let us take 
the differentials of the quadric with regard to #,,,...,, then the resulting 
equations involve none of the ca, but only the dn They reproduce in fact (xxxiv.). 
But the values of un... Y, found from (xxxiv.) are, we have seen, identical 
with the values of (xxxu.). Thus we have 











Pe en ~| R+ Dina Zait 4+ Rat Daga Ti+ h. + St R(q+ Ian Hai h E 
G SE P S 2 Se ? 
cc R@t Dyan 91 Hioilkuaun ge R@+ Dye da 1 
E l RGA gto Go E Latas Tyra E 1 Rg +2 Firo h 
SE R(g+2)tzq 01 : Rat 2ta Ta 8 l Riotiluaunn Oy de 


, BO ia C RAO: che SC 
EE EE Ke Si h, -+ R Ser e ho Së ER + R Së 5, hy ues ‘ ; f (XXXvi.). 
These give the most general form of a theorem proved for a particular case im 
‘Phil. Trans., A, vol. 187, p. 300, e (11.). Lf systems of q organs be selected with any 
arbitrary variations and correlations vut of complexes of n organs, then the mean sizes 
of the remaining n — q organs have precisely the same values as if the selection of all 
the systems of q organs had been to one size and not varied about mean values. The 
arbitrary variations of the selected systems about these mean values, as well as the 
arbitrary correlations, have no influence on the mean changes of the n — q organs. 
Returning to equation (xxxv.) we know that if the determinant 


| teen — — 
A — | Cad C1 Cy 9 . H . C1 gs C1 g+] e e e C1 in 
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| 
| 
| are. Kies es SEA KE) 
SE Geer, Lë, ea a Si Geen EC e (xxxvii ) 


Carito Care 941,38 + + © Cory Cottage ee © Corin 
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D) e HM e H e é e D e DH D e e $ H H 


Cn, ls Ca, 23 Ca, 3 S À $ Cn, KÉ Ca, +1 S $ $ Can 


be formed, 1ts constituents and not the linear terms in the exponential of (xxxv.) 
determine all the standard deviations and correlations. Let 3, be the variation 
after selection of the u™ organ; then if u be one of the selected organs $, = s,, if u 
be for one of the unselected organs 3, has still to be found. Let ta be the correlation 
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coefficient of the w and vo organs after selection, then 1fu and v be selected organs, 
Cu = Pw; 1 they be not both selected organs, then 1t has to be found. Let D,, be 
the minor corresponding to the constituent c,, in the above determinant, then by 
(x.) and (xiv.) 
: bes oe 
DS ae E (xxxviii. ), 
and 


Al, Ll E E 


Now let us write for brevity a, = D,,,/A and oa = D,,/A. Clearly as long as u and v 
are less than q + 1, o and #,, will both be known, t.e., are equal to s, and s,8,p,,. 
Now by a well-known property of determinants 


Cy, 41% 11 Co, 941% 19 Gaang T > Cg ga Sig Cpt 91% 91 be E Ong yin = 9 


Ci goði E Co, gotyo F C3, 94913 F > + + E Cog 49019 E Cpe igro%g4a Pe F Crge = H 


Cy, nO F Co 112 T C3 113 E Com Ry 9 F Cao 41,0], 9+2 aie Se Cm 0 a (x1.). 


Comparing these equations for the n — q unknowns Gun, 749.» + Gu With (xxxiv.) 
for finding gpu @ 42+. 27 u we see that they are absolutely identical if we change 
hi, hy... Agin the latter into ou, a)... Gu, Accordingly the solution is given by 
(xxxvi.), or we have 
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Provided v be < q + 1 we might equally well have used ou, %.,... Ga Hence 
we conclude that if v be < q + 1 and u be >q, 
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Or, substituting the known values of a, Gu... Gau, we have forv<q+ 1 and 


“>q 
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R (20), 8 


KR@) 8 7 Hoh 
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The next stage in our work is to find a,, when both v and u are > q, and also a,,, 
when v = u and both are > q. This is done at once by substituting the minors æn, 
Ayo) + + + by. +. Ay M the equations formed from the last p—q lines of the determinant 
(xxxvil. ). 

We obtain the following system : 


Cy 41 buy + Co g+1%A 02 SR C3 at1ffau + E, + Ey 41% ng + Carr +u, ot) + o + Cn, +12un Ges 0 
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An] + Con A St C3 Az Eg Se S + Cg, uA yy Ss Gott. A oi? + Ne, E Co uun Ss 1 


Cy nein AR Con Ana SG Czn n EE + Cy nny SG Gott, nn gH) T age F Cn, n Aun SC 0. 


These are identical with equations (xl), except that the equation with the 
coefficients Cim Con, - ++ Gan has unity instead of zero on the right-hand side. Hence 
we see that «n (v > q) will be the same function of am, Ay, . . . %. that æ’, is of h, 
ha... hy in equation (xxxv1.), but it will add to this a function of the last p—q 
system of c's, t.e., the Us 


Carlo Car gro + + + Cugei 


Datt Sgtigto e + + Cangt 
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Whatever this function may be we will represent it for the time by y,,; we notice 
that 1t is independent of the selected variations e, Sy, . . . Sp the selected means h,, 
h»... hp and the selected correlation coefficients pj, Ben... Daun It depends 
only on the characters before selection. 

We thus have 


ingen eae Et Drei 
SE) H ; Ye D" S qu 
A a Ye m ee C Ee E ek 
lit Vu E (Cl, o ul + R Cl T, u2 + + P (Do 7 
Now the system 0,1, Ar,» . . Gan can be found from (xli), since yw = Gan, whatever 


u and u” he. 
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Hence we have 
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Cy Op 
Here $, denotes a sum from p = 1 to p = 9. and H. a sum for every pair of values 
] } q, 2 
of p and p* out OF E EE 
When u = v we have simply 
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[t only remains to determine y,, and a, This we can do by putting all the ss 
zero, or selecting our q organs of one size only. We see at once that y,, and y. are 
the values of æ, and aw, that is, of Stu and SZ, when we select g organs of 
definite values and seek the correlation and the variabilities of two others, the a" and 
the oi. These values have already been found on p. 10. Or: 

Sli — Vue Se oimh (260). R(Q), | 

2 = Y =o k (ul Ro). a ee EE 
E E > | 
O Ge Sy R (U0 ki R (9), 


The notation of that page has been changed so that R (we) now stands for the 
determinant 
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R(uv),, is the minor corresponding to the constituent rw; B (uv), the minor 
corresponding to the constituent at the meet of the w™ column and u row; and 
R (q) the determinant with the last two rows and two columns struck out, 


For example : 
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We can now collect our complete results. 
The variability of a non-selected organ v >q is after the selection of y organs 
given by 


i R (uv), 
EEGEN 


E (lo > ëch? R (lo R (la SpSmu | | 
Lie E E SE ES 
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The correlation of two non-selected organs v and u both > q is after the selection 
of q organs given by 


ees R (a Etna Ro. se 
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The correlation between a non-selected organ u > q and a selected organ v<qg+ 1 


is given by (xl1,)% or, 


R (uhin Prv o) 


R (u) y SEN .. 
See = SO y S, Ee be l A ` a : (xlvii.). 


* The expanded values of these determinants are given, * Phil. Trans.,’ A, vol. 187, p. 294, 
NOT, CC,——A, Dp 
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Here in $5, : p takes every value from 1 to q, and in S,: p’ and p” every possible pair 
of values from 1 to q. Equations (xlv.)-¿xlvi.) fully determine all the required 
quantities and form the full solution of the problem of selection. Before we see the 
remarkably simple forms they take for simpler cases, we may draw some general 
conclusions of a most important character. 

In the first place we must distinguish between directly selected and what we 
have termed non-selected organs. It would be better to term the latter indirectly 
selected organs. Suppose the recruiting sergeant were to pay attention only to 
stature and seek to form a regiment of men of about 5 feet 10 inches. He might 
have a real range of stature about 6 or 7 inches, but he would strive to get men 
of about this height from the population. We will suppose that he did not consider 
chest-breadth, head-length, foot-length, lungs or any other character. The distribu- 
tion of these “ non-selected” characters in the regiment would not be the same as in 
the general population. Their means would have changed by (xxxvi.) and their 
variabilities and correlations be given by (xlv.)—(xlvu.). In other words, an indirect 
selection would have taken place. A selection by stature would change foot-length 
and head-length and indeed every other correlated organ. Much the same result must 
occur in natural selection. Ifit be advantageous for a species to have a certain group 
of its organs of definite size, falling within a definite range, and related to each other 
in a definite manner, then these changes cannot take place without modifying not only 
the size, but the variability and correlation of all the other organs correlated with these, 
although these organs themselves be not directly selected. Practically this means 
all the other organs, for so far one can hardly say with certainty that we have come 
across any two characters in an organism which are uncorrelated. Many of those 
investigated are highly correlated, all appear to have some correlation, even if 14 he 
very small or negative. 

We may therefore conclude as follows :—- 

(a) The selection of any complex of characters or organs in an organism changes 
all the other characters and organs not directly selected. 

(0) If the change in the complex be continuous and progressive, the other 
characters will continue to be modified until the change in them is so considerable 
that selection begins to act directly upon them also. 

(c) The changes noted here are not confined to the average value of a non- 
directly selected character and to its variability ; the correlations between non-directly 
selected characters and the correlations between directly and non-directly selected 
characters are also both changed. 

(d.) If local races have been produced by selection from a common stock, it will 
be impossible to look upon correlation as a criterion for species. Every selection will 
modify such correlation, and it has no greater fixity than either type value (mean) 
or variability (standard deviation). 

The whole of these statements will become more manifest as we apply our general 
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theorem to special cases, but we must note that if two organs were uncorrelated with 
each other, it still might be possible by selecting a third, or a third and fourth, to 
produce correlation between them. Further, by selection of one or more organs, two 
non-directly selected organs can have their existing correlation increased, lessened or 
even changed in sign. 

(4.) A primary difficulty will of course arise in the case of natural selection. How 
are we to determine which are the directly and which are the indirectly selected 
organs? With artificial selection by man, we know which organs have been selected 
fairly well; attention has been paid to colour, size, proportion of parts, Ae, Even in 
the case of the medical examination of the recruit, 1t is chest, lungs, heart, stature, 
&c., which form the basis of the acceptance or rejection. If the head or foot be not 
absolutely deformed, little if any attention is paid to them, so with hair-colour, 
probably eye-colour, and a mass of other details. No doubt the direct medical 
selection indirectly selects these, but we could roughly class the selected and non- 
selected organs or characters and investigate the changes in the correlations of the 
latter owing to the indirect selection. But how are we to form these classes in the 
case of natural selection ? 

The investigations may look difficult, and even from the standpoint of arithmetic 
appalling, but it seems to me that the differentiation of organs into directly and 
indirectly selected classes ig the keynote to the problem of evolution by natural 
selection. 

Let us look at a simple case and see whether it will throw any light on the 
problem of distinguishing between directly and indirectly selected organs. Suppose we 
have two organs only, with means m,, mo, standard deviations ex, oa, correlation 7, 
and let the first be selected so as to have a mean value m, + Ap and standard 
deviation s, Let % be the standard deviation of the second organ and fti the 
correlation of the two organs after selection, and m, + x’, the mean of the non- 


selected organ. 
Then by (xxxu.) : 
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In other words the regression coefficient of the nou-selected organ on the selected 
remains unchanged, while that of the selected organ on the non-selected will, as a 
rule, be widely modified. | 

Further, let X, be the mean value of the second organ before selection corresponding 
to a value H, of the first; let M, and M, be the means of the organs after selection, 
and Y, be the mean value of the second organ corresponding to a value K, of the 
first. Then the equation to the regression line before selection 1s 


O; - O. 
BC Z Pig SS H, Se Ho, do OR SC Vey, 
R + o SE S o 


and after selection it is 
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But this 1s identically the same line as the regression line before selection. Hence 
not only the slope (regression coefficient) of the line, but its position is identical, and 
we have the following result :— 

If two local races have been evolved from a singie stock by the selection in different 
ways of one organ only, then the regression lines for the two races of any non-directly 
selected organ on the directly selected organ will be the same in direction and position; 
but the regression lines of the selected organ on any non-selected organ will differ for 
the two races.* | | 

Of course the means, standard deviations and correlations, not only of the selected 
organ but of all the non-selected organs also, will probably have changed. It is only 
certain of the regression lines which remain unchanged and serve as a criterion to 
enable us to distinguish between directly and non-directly selected organs. 

Of course the problem in Nature will not be as simple as this, for differentiation of 
the two local races may have arisen from the selection of more than one organ, or 
may have arisen from the selection of two different organs, but the illustration will, 
I think, indicate the nature of the investigation we are proposing. 

We can easily generalise our theorem by considering the form of the selection 
surface given on p. 12. Any result obtained from (xxxv.) which does not involve any 
of the ës will be a result unaffected by the selection that has gone on. Now to 
obtain a regression equation we put any number of the op equal to constants, to h's 


* The geometrical interpretation in this simple case that the regression line is unchanged is quite 
obvious, and, indeed, may serve as a proof. ` 
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say, and find the “centre” of the quadric of the remaining og, the co-ordinates of 
this centre, expressed in terms of the Aa, are the regression equations. Now it will be 
clear, that if we put all the selected a% equal to Ae, the differentials of the quadric 
with regard to the remaining or non-selected ais can contain no da or the coefficients 
of the regression equations thus found will not be modified by selection. 

Further, we might have given not only the selected organs, but any number of the 
non-selected organs constant values, and the resulting regression equations would 
involve only the ds and not the cs. 

Hence we have the following general theorems : 

(1) Lf an organ has been modified only by indirect selection, then uts partial 
regression coefficients on any complex of other organs, however large or small, provided 
it includes all the directly selected organs, will remain unchanged by the selection. 

(D. The same organ in two different local races which have been derived from a 
common stock by the selection of two complexes of organs, some of which may or may 
not be common, will, if tt has not been directly selected in either case, give the same 
partial regression coefficients for any group of organs which includes the members of 
both complexes and any number of non-directly selected organs besides. | 

1f the partial regression equations have changed coefficients, then we cannot at 
once determine whether— 

(a) We are dealing with a non-directly selected organ, and have not included all 
the directly selected organs in the group upon which we are calculating the regres- 
sion ` Or 

(b) We are really dealing with a directly selected organ. In this case, we have 
also certainly not included at least one directly selected organ in the regression 
group. 

Theoretically, however, (1.) and (i1.) suffice to find out which, if any, are the non- 
directly selected organs in the differentiation of local races. Practically, however, 
the number of organs and characters may be so great, and our ignorance of those 
probably selected so complete, that the arithmetic of determining so extensive a 
series of partial regression coefficients may be quite beyond our powers. Still, where 
the divergence between local races is not too great, and the source of the differentia- 
tion not too obscure, it is probable that the above theorems will lead to results of 
great interest. ™ 

Without laying too great weight on these theorems, I would still venture to 
suggest that if the criterion of a species be the discovery of any numerical constant 


* Mr. L. BRAMLEY-MOORE has been working with this end in view at the long-bones in man. But 
even here the direct selection of parts of the vertebral column—for which, at present, we have no correla- 
tion values either among themselves or with the long-bones—and of the hand and foot, which 
Dr. W. R. MACDONELE. has just shown, are very highly correlated with the long-bones, may render 
nugatory all attempt to ascertain which, if any, long-bone has been only indirectly, or, at any rate, least 
directly selected. 
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or group of constants, which is the same for all local races, then these constants must 
not be sought in the values of mean characters, degrees of variability, or of correlation, 
but in a system of partial regression coefficients, and the discovery of these is therefore 
of first class biological importance ; it is the classification of the characters into directly 
and non-directly selected groups, %.e., it 1s the discovery of the modus operandi of the 
factors by which the differentiation has taken place. We are a long way from 
solution yet, but we may venture, perhaps, to admit a fait glimmer of light in the 
direction of what might seem the culminating problem of the mathematical method 
as applied to evolution—-the piecing together by quantitative analysis of the stages of 
descent. 

(5.) I will take now the application of the above results to simple cases ; but for the 
benefit of those who cannot easily follow the main principles of our investigation 
through the stages of determinant analysis, I will prove directly the proposition that : 
ihe selection of an organ A alters the mean and variability of a correlated organ B, 
and also the correlation between A and B, 

Let the correlation surface for the two organs be 

N di ES AU 2 Zait A te ) 


o H AAA ie e Nal) oad) o Ar) 


REN H Ee SE we 
270,0, Y 1 TE Zus 


D 
d 


where N is the number of individuals in the general population before selection, and 
the subscripts 1 and 2 refer to the organs A and B respectively. 
Let the distribution of the population after selection of the A organ be 
(e, — hy 
H EL 
¿== ae Y, 
where N — n is the total destruction, 4, the mean and s, the variability of the 
population with regard to A after selection. Before selection this distribution was 
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Hence, the selection being random with regard to the array of B's corresponding to 
any A, we have for the surface after selection 
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for each array must be altered in the ratio of the corresponding 2” to zy 
This gives for the surface in full 
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Integrate this for every possible value of x, and we shall have the distribution h, and 
s, of x, or B after the selection of A, After some reductions we find for the frequency 


(dx, 


where 


KEN Ss ge (1 ar y KR | a) fer è . ‘ > . ‘ 3 (xlix.), 


which gives the standard deviation of the indirectly selected organ, and the “ centre ” 


of this organ is given by 


hs — GES h, . » e a e . > 2 . 2 . (1.). 


di 


These are special cases of our results (xlv.) and (xxxv1.) above respectively. 
Further, returning to the correlation surface (xlv111.), the coefficient of correlation 
tp is the coefficient of xx, divided by the product of the coefficients of x? and g. 


Hence we find 
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Let zua = sin Bus, Dua = SIN due, then we have 
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This shows us that du, decreases with a, or that Ca decreases with s,, that is to say, 
the more intense the selection the less is the correlation. This in broad terms 
demonstrates the general principle that intensity of selection connotes a lessening of 
correlation. Tt is this principle which very possibly accounts for the fact that the 
more civilized races of man appear to be not only more variable but more highly 
correlated than the less civilized, among whom the struggle for existence is more 
intense. It may, perhaps, also account for the skeletons of women of the civilized 
races having their parts more highly correlated than the parts of those of men.* 
Lastly, it may well throw some light on the markedly plastic character of races 
which have been stringently selected with regard to one or a few organs only. 

As an illustration, suppose that the correlation between femur and tibia were -7, 
and let us investigate what would be the effect of reducing the variability of the tibia 
by direct selection 50 per cent. We find at once on substituting sie: = '5 and 
ra = 7 in (h.) above that ty = (44, or a reduction of about 37 per cent. This will, 
perhaps, be sufficient to indicate what immense changes must be made in the 
correlation of highly correlated organs whenever selection, artificial or natural, is 
stringent. It is important to notice that the change in the size of the organ in no 


* See Lux and PEARSON, * Roy. Soc, Proc., vol. 61, p. 354; and Les, “Phil, Trans.,’ A, vol. 196, p. 231. 
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way influences the change in the correlation between organs, if the distribution be 
normal,* the change depends only on the stringency of the selection. Breeders who 
select. by the size of an organ only are in that case very likely to reduce the variability 
of the organ in the selected group by far more than 50 per cent. Accordingly, 1t is 
not to be wondered at if they, to a great extent, destroy the correlation between the 
selected organ and other organs. This destruction would appear as a want of relation- 
ship between parts, possibly as a want of fixity in type. 

By means of equation (Un. ) t} can easily be found from 7), and p=s,/o, All we need 
is a table of trigonometrical functions. We observe that r, is always of the same 
sign and less than tj, For many biological purposes the following graphical con- 
struction gives quite sufficiently accurate results. Let CAB be a quadrant, say of 
10 centims. radius, and take the point P on this quadrant distant PN = 107, 


from AB. Take QN zz: PN, and let AQ meet the quadrant in R, then RM the 
P Š 
distance of R from AB = 10r,,, and consequently determines t, If the figure be 
drawn on decimal paper the determination of t} 18 peculiarly easy. 
Graphical method of finding correlation between organs A and B after selection has acted on A, 
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— { y P 
In the above example au = PN = “7; syo = =- = 4, and therefore QN = n/p = :35. 
pP 


te = Ry = RM = :44, as before. See p. 23. 








* Tt will do so if the distribution be skew, see ‘Phil. Trans.’ A, vol, 191, p. 231, 
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Further, taking AP as our unit, AN? = 1 —7,2?, and QN? = rg 


from (1x1x.) : 


Hence, 
3? = 07 (AN? + QN?) = 0AQ?. 


Therefore if Ao, in the diagram be taken equal to o, and QS, be drawn parallel 
to Ro,, we shall have A, = Ze, or we can scale off the reduced variability. 

Thus the diagram enables us to see at a glance the reduction in correlation and 
variability. 

(6.) Let us now write down the results when an organ A is selected out of a group 
of three organs, A, B, C, whose constants are marked by the subscripts 1, 2, 3, 
respectively. Let Ay = gue, and be represented, when required, by cos zw, Then 
we find from (xlv.) (xlvii.) : 


, Ee Og } ELA Og eee 
Xx DAR Tig Nee Gë Lg = 1) 13 h, e e e e D e e (L111. ), 
Se o 
E d 29073 
Za = 0931 — (1 — =) mé] = os {1 — sin? yı cos? Oig}? 
i D GE 
1) L 


Se ña 8 ae 9 z e 1 
l; 








E Haig _ COS wu COS Oj 
| {1 a E — sin? Xı COs” E (lv ] 
E EN in, 
ti = j~ 





E A (1 — Du 2) r} Ya E sin? X1 cos? Cp 


9 9 
S 8 
S 3 $ 
(Paz — Pig" (1 = A REES 


97 oO} 


CH Vlad dal EY ia, 


COS Oa; + sin’ x; Cos Oi Cos Oig 


AN (1 — sin? y Ae cos? b) (1 — sim? yı COS 18? d. 











where, as before, we write 7,, = cos 6,,.. Let us also write t,, = cos @,,, and 


SIN wi COS Ha = COS Ayo, SIN vu COS 0), = COS (ty. 


Then we can replace the above results by 


29 == Fy SID 0, 23 = os SIN dhs, 
COS ©), = cot y, Cot Ayo, COS Oj = cot x; Cot (ys, 
COS Ay, — COS yy CUS UU E 
COS Dos = : E aw ` y CS D D D D ` D D (Lon) 
i SIN tja SIN Ajg 
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These equations admit of easy interpretation by spherical geometry. 





Let P be the pole of the great circle DEFG. Take DG = bap, DE = 0,,, GF = fu 
Join P to E and F; let the small circle of radius y, round P meet PD, PE, PF, 
and PG in d, e, f, g respectively. Draw the arcs De and Gf Let the small circles 
with centres D and Gand radu De and Gf respectively meet in Q. Join DQ and GQ. 
Then the quantities required are : 


Si sim DG), So = su GO). 
D = cos Dek, == cls 16. lag = cos DQG. 
es S T TT 
For DE = 0, we = AEN Dike = oe 
hence : cos De == cos 0, cos He = sin y, cos 0,4 == cos gu, or a, = De; 
similarly Ons Gf. 


Next, cos DeE=cot De tan eE =cot os cot yp or 2, DeE=0,,; similarly < F/G=0),. 
Lastly, from the triangle DQG : DC) 0 5, QU and DG bs. but 


cos DG = cos DQ cos QG + sin DQ sin QG cos DQG, 








en cos DQG == ~ SE = COS Oy; or L DQG = Oy. 

Thus all the relations can be expressed in terms of the sides and angles of a simple 
system of spherical triangles. For the degree of accuracy generally possible in 
biological and sociological investigations these triangles can be solved by a spherical 
trigonometer, such as that sold by Krurpn, of Prague.“ The changes, however, 
which tẹ, undergoes for various values of ria Zus, Tog are, indeed, far more difficult 
to appreciate as a whole than those of t}, or rı In order that they may be followed 
easily, and in order to solve directly to a degree of approximation sufficient for many 
practical purposes problems in the influence of selection on correlation, my assistant, 
De. L. N. G. Finox, has kindly drawn up the tables which accompany this memoir. 


* Tt will suffice fairly well for all but a few special values of 71», roa, F31 
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Of course to bring them within any reasonable compass we have had to limit the 
values taken. In the first place we have considered only eleven grades of selective 
stringency given by 


s/o, = 0, 1/10, 2/10, 3/10, 4/10, 5/10, 6/10, 7/10, 8/10, 9/10, 1, 
the corresponding values of tẹ, in the tables are entered as 
Eo. E Bae Rese Ro Re Re, Ro Eo Rig 


The tables are calculated for 7, both positive and negative, but 7, and rig are 
always supposed positive. If ra and ri be both negative, then ta will be the same 
as if they were both positive. Tf ra and 7,3 be of opposite signs, then all we have 
to do is to look out te in the table in which 7, has a sign the reverse of its actual 
value, and having found the corresponding value of ta, then change its sign to obtain 
the actual coefficient of correlation after selection. This follows, if 7,, be the 
negative coefficient, by writing : 


4 4 4 i 2 
Taia e ON Yie 


T E EE 
TT {lne siny {1 — yy sin? y ` 








T V ne sint) i — n E aF 








Lastly, it would clearly be very laborious to tabulate t,z for a very great series of 
values of ria 713, Tan, Accordingly a selection had to be made of these coefficients of 
correlation. They were given the values 0, :25, '5, 28. and 1. These may be spoken 
of as zero, small, medium, large, and perfect correlations, and the ranges 0 to '25, 
:25 to ‘5, '5 to ‘75, and 75 to 1, as the ranges of little, moderate, considerable, and 
high correlation respectively. There would thus appear to be 15 combinations of values 
for rig, Ta: these are given in the key to the tables as (a), (b), (c), (dl)... (m), (n), (p), 
see p. 63. If these 15 values had to be combined with the 10 values (5 positive 
and 5 negative) of 75 and the 11 values of s,/o,, we should have 1650 entries in our 
tables. But this number is much reduced by the consideration that the expression 
L — Py — fu — Mig” + 29113149 has for the real correlation of three characters to be 
always positive. Yo can also never be greater than unity. Accordingly all values of 
Fog, igs 719, Which do not satisfy these conditions, have been excluded from the tables ; 
they cannot arise in nature. A few impossible values of tẹ have been included in 
the tables, but these are placed there solely for the purpose of finding by interpolation 
values of Ga, which are less than unity. The following purely hypothetical illustrations 
of formula (lv1.) and the tables will serve to indicate their use. 

Illustration I.—Suppose the correlation of tibia and femur with each other to be 
:8, and of both with the stature to be ‘6. How would their correlation be altered if 
the variation in stature were reduced by selection to half its present value ? 

Let s,/o, = p as before, and suppose taz = R; then let p, fio, 7,3, Pag be the values 
of the constants next below the required values occurring in the tables, and giving 

E 2 
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{oy == R; let R’ be the true value of 1.3, corresponding to the values q, + òu, tig + Oro, 
Tis + OP 13, Tan + 873. ‘Thus we have, as far as first differences : 


R’ = R + 10 (A,R) 84, + 4 {(A,,B) ër + (A, R) 875, + (4, R) ër), 


I case == Ms E 00, POR == ol, == Ue 
Gin = ‘05. Further, we look up Table IV. (a), and the nearest case is (7) under R,, 
which gives R = 67105. We then see that (A, R) + (A, R) = difference between 
(7) and (m) cases = — ‘21455; (AR) = difference between result in IV. (a) and 
V. (a) = 82895; and lastly (A, R) = difference between R; and R, columns of (7) 
row of Table IV. (a) = ‘00535. Thus we find: 


P13 


R’ = 671050 + :002675 — 04291 + ‘06578 = '6966. 


The value by straightforward calculation of formula (lvi.) is ‘6981, the two results 
giving substantially the same value ‘7. Thus we see that such a selection would 
reduce the correlation of tibia and femur by 12:5 per cent. 

Illustration (I1.).—Suppose the correlation of humerus and femur to be ‘5, and of 
those with stature to ‘7 and :8 respectively. How would the correlation of humerus 
and femur be modified by a selection of stature given by s¡fo, = "5? 

ing tie Case. eye Oe O Ee 
Bi = 0. We turn to Table III. (a) and take out (el under R,, which gives us 
R= 38192. We have A „R = — ‘1841 and A, R= — :04185, whence we find 
R’ = :1636, but the differences of the table are too great at this point for the 
result to be very trustworthy.” Suppose we take y, and 7, as before, but r = "75, 
r = 75, and therefore R, to be found from (m}, = 1351; then dr, = — 05, 
ôr = "05 and A, R= — 1841 as before, A R = — 2995. Hence we deduce 
R’ = 1120. The mean of these two values of R’ is ‘1377, and the true value 
calculated from (lv1.) is R’ = +1395. Taking '14 for the practical value, we see that 
the correlation of humerus and femur has been reduced by this comparatively 
moderate selection of stature upwards of 70 per cent. ! 

Illustration (111.).—Suppose a case in which humerus and femur were not correlated, 
but that both were correlated ‘7 with stature. What would be the effect of 


selecting stature with the same intensity, Ge, 8/0, = "A? 

Our best results from the tables will be to take R; (m) from Table I. (a), which 
gives R = — ‘7297. We have then zs = rig = ‘75, hence ëss — 05, òr = —'05, 
An R is to be found from (m) and (4) and = — ‘3193 = A, Rh, and 

R = — ‘7297 + (A, R) ôr + (A,, R) äis, 
= — 7297 + :1277 = — :6020, 


The actual value by formula is — :5810. 


* Second differences ought to be used, and the process indicated is practically equivalent to using 
them, 
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Now this again is a remarkable result; by selecting an organ correlated with 
two others, neither of which are correlated with each other, we have produced a 
considerable correlation, and what is more, one of a negative sign. 

In other words, 1f humerus and femur were unrelated to each other, but were 
related to stature, then a selection of stature would result in men of long femur 
having a short humerus, and vice versd. | 

Illustration IV.—Suppose the correlation between greatest length and breadth 
of the skull to be ‘25, and between both and the auricular height to be ‘5. Now let a 


1 


stringent selection, A = 315, of height take place. What modification will there 
S F 0 H p l 


1 


be of the length and breadth correlation ? 
Table IT. (a), R,, case (7) gives us at once the result— 


Los — "0038. 


In other words, the correlation between length and breadth would be sensibly 
destroyed by such a selection. Thus correlation can be created or destroyed or 
reversed by selection. 

The above illustrations, hypothetical though they may be, will suffice to indicate 
how entirely dependent correlation is upon selection. We must look upon coefficients 
of correlation, in fact, as just as much the outcome of selection as coefficients of 
variation, standard-deviations, or even the mean size of organs. No selection can 
take place, in the sense in which it has usually been understood to take place—+.e., 
by a change of mean and of variability, without at the same time the means, 
variabilities and the correlations of all correlated, but not directly selected, organs 
being varied. This is true whether the non-selected organs be initially correlated or 
not among themselves. We must always bear in mind this all-important funda- 
mental conception, that natural or artificial selection, or even random sampling, are 
in themselves active factors in the modification (e, creation, destruction, or reversal) 
of correlation. Thus not only is the impossibility of the constancy of correlation 
for local races obvious, but the primary importance of insuring that our samples are 
representative, and not accidentally selected samples, in all observations or experiments 
on heredity, homotyposis, or organic correlation becomes more and more manifest. 
We must not lay too much stress on two heredity constants— differing, for example, by 
more than the probable error of their difference—unless we are convinced, which 
practically it will be difficult to be, that all modification of correlation by unin- 
tentional and unmarked selection has really been avoided. 

(7.) Let us now take the next most simple case. Jf A, B, C, D be four mutually 
correlated organs (in either the same or different individuals), and a selection take 
place of A and B, to find the changes in the characters of the non-selected 
organs. 


o 
e) 
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Let subscripts 1, 2, 3, 4 mark the organs A, B, C, D respectively; let Au, ha, 
gu, ën, Pig be the constants which determine the selection of A and B; and let us 
apply (xlv.) to (xlvi.). Here the coefficients will be obtained from the set for q = 2 
in (xxxiii. ) and (xliv.). Hence we find : 

















? D a a a 9, 9 l ga wees MP H 9 
O K — Pig — Mig — Pog + Zëesfag ete Tig Coste (5 ae Pog — Yigg ch 
ein = Oe 1 9 1 2 1 o, 2 
L - — h SE \ oy "Tu Ta] 
Pig Tas Tog — Deia $189 bt 
+ 2P13 > GER EE VIL. 
Li — P — Poe + Hauf as e ET A re 
Sti ee 12 14. 24 197 14794, Pg — Poata D E Tii Ki 
i See A EECH PRS o 
S AC ER a? O) L= ng da 


— Vea" 1a Po — Tiatia $159, 


le 





o 
pa) 


12 


8 
+ 2P12 El 


pe 


4j@ 9 a e] 4 fy? kl a a? a 3 
Pa (d — Pig) — ul — Toston + Pra (Pisa + 114723) 


A 
a E ee 








Zadata, = 037, 


l re 


f 
































4 Pigs — Tagtig Pig — Tatra [$ V E Vaz — Vigia fa — Dm (SV 
a a a ER D 2 q D a 2 mm ft MP ` 
Tp rig — Togia Fog Tu ES Pra Tag Tog — T1129 \ $189 (Ix ) 
Pu POS Ee 8 : 
13 2° 93 1 23 1915 2 
gata = HC eae + = SE i E ME 
143413 173 l E o Den. Pia 5, (1x1.) 
PR TES Po E $ o 
Sia, es My Ga %1 Pos 12 9 Ën 
D i A 2 O. e e e o ie e 
14414 1" 4 eS m2 o Su 1 A 9 12 o, ( ) 
; Tox — Ty27"13 82 T13 — T1993 S1 see 
Sos = $203 KE 2 Ss -+ KE EE 19 E ° > 0. (1x111.). 
12 a 12 1 
d a Pos — Tel 14 Se Pri — Tigro $] os 
Si. = Shey l n P -+ qu a E (Ix1v.). 
— 48 2 Ce E 1 
Finally, for the change of means of the non-directly selected organs, we have : 
Pro — PioPos E Pos — E Ta 
fee as E S 12°23 3 23 Sa os AS E 
a BE ES WË h 4- HE hy s o e (Ixv.). 
Lg a L— P 0 
r Volo, O Poy — Mit, O : 
l4 1994 Ty "oa Tu 1 
al, = HSH ZA BI. (xvi). 
> (md 0 Ee Li, a . l 


If we write (lxv.) and (lxv1.) in the form— 
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v= = By ` SE + Baz a 


E , (lxvii), 
wy = Br. T Ba Se? 


the Gs are the partial regression coefficients, aud the whole solution can be expressed 
in terms of them. Thus : 


5189 


Zell = gell d — Bisis — Baatz + Br (2) + Bay” (2) Se pff ` SC Di (Ixviil.). 


y Oy 


Zeta = 030) d 34 77 Bis" — Bastos + ME 3] + Ba E a] 
+ pes (Biba + Bai, a el ee, (xx), 


82113 = — SIS {Bis = += + P1903 a S ‘ s . d R (Ixx.). 


2 


Thus the whole series of results can be easily calculated, if the regression coefficients 
are first calculated. 

I may make some remarks upon these results. A formula equivalent to (lxviii) 
was first given by me in my memoir on “ Heredity, Panmixia, and Regression” 
(‘ Phil. Trans.’ A, vol. 187, p. 303), and used for certain problems of inheritance, and 
conclusions drawn from (lxix.) or (lxx.) have been cited or- indicated in other 
memoirs. 

Some interesting results follow at once. If the selection be very stringent, s,/o 
and s,/o = 0 sensibly, then all correlation between a selected and non-selected organ 
is destroyed. But 





A "aal = J LS EE EE H Pi (7 Pigo E figa 3) a 7 (Ixxi.) 
a vU = r — Pag — 743" + 2191130 95) tt Ty gy ne an Zait) ` 


This is what I have termed a partial correlation coefticient—1.e., the correlation 
between C and D when fixed values are given to A and B. So far as I am aware, 
such coefficients were first directly used by Mr. G. U. Yure in certain economic 
problems.* They are of very considerable interest, but for natural or artificial 
selection are not quite so important as the generalised form (lxix.), for we generally 
select about a mean value, and not absolutely at it. 

It will be noticed that the coefficient of correlation of two non-selected organs 
differs from the corresponding partial correlation coefficient by terms of the square 
order in s/o, but the coefficient of correlation of a selected and non-selected organ 


* «Roy. Soc. Proc.,’ vol. 60, pp. 485, 488; * Economic Journal, December, 1895, and December, 1896. 
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differs from zero by terms of the first order in s/o. Hence, when selection is intense 
(s/o small), we may neglect the former as compared with the latter, and we have 
thus the basis of a method of approximation very useful 11 some cases. 

I have not yet succeeded in giving a geometrical interpretation to the above 
formule, but have little doubt that it would be quite easy if the “ spherical” 
geometry of four dimensioned space were more familiar to me. It will suffice to 
observe that it is easy to find cases in which the correlation of a non-directly 
selected organ with a directly selected organ, or with another of its own class, can 
be reduced, destroyed, increased, or reversed. In fact, all our previous warnings as 
to the caution necessary in avoiding unintentional selection in collecting material for 
testing correlation remain in force, and, indeed, are emphasised. 

The following illustrations will indicate the kind of problems which may be 
attacked by such formule as the above :— 


(8). ILLUSTRATION I.—Influence of a Selection of the Long Bones of the Leg on the 
Long Bones of the Arm, and on the Relation of the Leg to the Arm. 


A numerical example will throw light on the application of the above formule, and 
effectively illustrate the manner in which a selection influences size, variation, and 
correlation. 

Consider the long bones femur, tibia, humerus, and radius, indicated by the 
subscripts 1, 2, 3, 4 respectively, and let m,, my, Mg Mm, be the mean values in 
centimetres. Then the following numerical values are given in a memoir by Miss 
ALICE LEE and myself :—* 


French ¿. Aino d, 
my 45:23 40°77 
Mə 36°81 33°89 
Ng 33°01 29:50 
Ma 24:39 21:55 
Ti 2°30 1:90 
O) 1:80 1:67 
T3 1:54 1:34 
Ta KAE 1:06 
112 -806 "827 
113 :842 :858 
Pia 744 "189 
Tog :860 «745 
Toi K 780 i 865 
134 "HAD 176 





* “Roy. Soc. Proc., vol. 61, p. 343 et sey. The correlations are not worked out for exactly the same 
lengths in the case of the two races, but the numbers will serve quite well for the purposes of illustration. 
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Now let us select from the French population a group having the same 
characteristics of the long bones of the leg as the Aino population, and then compare 
the characteristics of the arm bones of this selected group with those of the Aino 
population. 

Our selection is given by: 


h, = 40°77-45°23 = — 4:46, 
h, = 33°89-36'81 = — 2°92, 
S100 pe = "827, 


S e 
The following constants must now be determined arithmetically : 


Sua = 190/237 S802, 
Si Onc OH 180 978. 





ae r Pos E EA 
e Ce SE 23 Se SE 
E — “426, 13 — ‘518, 





ES zz 
E 2 au 329, Poa Proa -515, 


A o e + 2719718193 





= 106, 


E E = 
Lom De + a 
oe 2 
Lie 


== :354, 


2 Ki 3 Aë > A nc a + A? kl + 3 
Poy (L — Py) — May" Tosfa + Tya (713 o4 + "sl 


ee |): 


If x’, and x’, be the mean humerus and radius ot a femur-tibia selection from the 
French population, we have from (lxv.)— 


e = 33:01 + :277h, + 4483h, 
al, = 24°39 + :162h, + :335h,, 


These would give the effect of selecting any femur and tibia defined by h, and A 


from the mean values of the humerus and radius. For the particular selection 
indicated above : 


Ve == 30°48, = 2269. 


both of which are about a centimetre in excess of the Aino population. By selecting, 
therefore, from the French, a population with a mean leg like the Aino, we should 


still find the average arm of this population some two centimetres greater in length 
VOL. CO,-—A. F 
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than the Aino. The variabilities X}, and 2, of humerus and radius for a population 
selected from the French by femur and tibia are obtained from (lx111.) and (lix.). 
We have: 


Y ~J 
¿ 1 sl 
S S 


S/o = 196 + 181 = + 268(-2) a22 py x 221, 


Fr) Oa ; O71 


Soe = 354 + 108 (= S sl + ‘265 dë y gt So TE 
aie | 3 2 


a, oT 
These give for the particular case : 
$= 1839, yell 


Turning to the correlation of humerus and radius, we have by A 


dE EE 140 (4 SE SR 267 | Sal + :389 pro Sé 


giving for the particular case : 


ta = 799. 


It will thus be seen that if we selected from the French a group with the same 
variabilities and correlation of femur and tibia as the Aino, the variabilities and 
correlation of the humerus and radius of this group would not be very different from 
those of the Aino. On the other hand, the correlations between upper and lower 
members would be very significantly different. 

Generally we have by (lxi.) for selection from the French : 

RER “Lt +518 p me 
Po, 


A ds 
Lia = e 
E 


Lu e laan? = a DES Paes = 


|! 


Tag 21 am + 426 po SH 


ty = E 2 + 829 pa 


Va GEI E 694, 
Yo, = = "HAN. Lan = 768. 


These are all smaller than the corresponding French values, the selection has 
reduced the correlation, but the Aino population has in all the cases but 723 a greater 
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correlation than the French. We must accordingly conclude that by a leg selection 
from the French aimed at reproducing the proportions of the Aino leg, we should not 
obtain an arm equivalent to the Aino arm. The divergences are indicated in the 


accompanying table :-— 





SE be E an A ua ae A Ee et 
NEE A E A er tere ett nt A 








Selection from the French. -© ÅAïno. o o 
| 

Mean humerus | 30:48 29:50 33:01 
EARL 8s chal oe 22269 21:55 24°39 
Variability of humerus. STEEN | 1:39 1:34 1°54 
S rau <2 e e A 1:11 | 1:06 Pary 

Correlation of humerus and radius . | "799 776 :845 

g , femur and humerus . | "819 ‘858 GER 

` , femur and radius . | "694 ‘789 "744 

We , tibia and humerus . | HAD 745 :860 

E „ tibia and radius ‘768 :865 780 





A A EE EE En 


There is, of course, no special reason for supposing that the French and Aino 
differ merely by an evolution which has acted by selection of femur and tibia. We 
might have obtained a race out of the French more nearly akin to the Aino by a 
selection of femur and humerus, but the process would numerically be exactly 
similar. The particular illustration here chosen is taken merely as an instance, to 
indicate how the methods developed in this memoir enable us to ascertain with 
quantitative certainty how far racial differences may be due to the more or less 
stringent selection of a limited number of organs in the one race. 

If we consider that local races have been differentiated from a parent stock by 
the selection of the chief or more markedly divergent organs, then we have in 
processes such as that just illustrated a method of ascertaining, at least tentatively, 
whether two races are to be considered as merely local varieties, and further the 
particular organs through selection of which the differentiation has taken place. 


ILLUSTRATION 11.-—Infuence of a Selection of Femur and Humerus in Modifying 
Stature. 


The following data have been calculated for me by Miss Arron LEE from RoLLETS 
measurements on the French :—* 


* They have been undertaking, with the view of determining more scientifically than appears to me 
yet to have been done, the mean stature of a race from a measurement of the long bones found in burial 
mounds, &c. ROLLET'S measurements are given in ‘De la Mensuration des Os longs des Membres,’ 
Lyons, 1889. I hope shortly to publish a memoir on the subject. [The memoir in question was published 
in ‘ Phil. Trans., A, vol. 192, pp. 169-244, 1898.] 
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3 e Y 
EE O eee g 
Mean stature,m . . . . . . . .! 166-26 cms. 154:02 cms. 
Variation of stature, 07 . . . . . +] DDO a, | Ee 
Mean lemur, Me "e EE 2. AR AË -, ARATO eee 
Variation of femur, og . . . . .. . | KEE 2720, 
Mean humerus. Mg . CR E EE, .%, SE T 
Variation of humerus, 73... . . «| od Jo | 103" 
Correlation, femur and humerus, 723 . | -842 "872 
ge stature and bumerus, 713 | ‘809 SE 
E stature and femur, riz | ‘Sil ‘805 
j 








Now let us select from among French males a group having the same variability, 
correlation, and mean size of humerus and femur as French females, and let us ask 
how this would alter the variability (3,), mean size (M,) of stature in French 
males, and also the correlation between stature and humerus (tj) and stature and 


femur (tis). 
We have at once from the second column— 


h, = — 3°66, ho == — 3:24, 


e 


KÉ == 2:26, Sx = 153, 


e 


whence we find, 





TRE BS a ‘447 Tyg = EE ZER 
1 — reel í | L = To ) 


From (lxv.) we deduce, 


M, = 166:26 + 1:037/, + 1:546hy, 


This formula gives the stature of any group of males selected from the French, 
and having their femur and humerus respectively hy and hz centims. longer than the 
average. 

For the special selection referred to above, hy = — 3°66 and hg = — 3:24, hence 


M, = 166:26 — 3°79 — 501 = 157:46, 


This example shows us that if we selected French men with the same femur and 
humerus as French women, it would be the selection of the humerus which would 
contribute mostly to the reduction of stature—a somewhat singular result. Further, 
such a selected group of French men would be still some 33 centims. taller than the 
average of French women (instead of about 124 centims.). Probably had we 
selected the tibia as well, the greater portion of this remaining advantage in height 


would have disappeared. 
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To find the variability in stature of the selected group we must use (lviii.). We 


deduce : 


9 A s \” e 
ZA o EZE 200 (* + 187 (2) + :387 23 pt. 


In our particular case this gives: 
$ =o, X 987 = 5:43. 


The actual variability in stature of French women is measured by a standard 
deviation of 5°45. Hence our selected group of men would be sensibly equally 
variable with French women, as far as absolute variation 1s concerned. 

Lastly, from (1xi.) : | 

GE 
EE [a47 E SH 


E 


EI ER 
Luz Ss S (um 8 + "447 = pas} 


Ox 


which give in our particular case : 


=> “8011 = “81, 


sec) 


ts =3 x 8017 = 812, 


Such a selection, therefore, would accordingly only increase insensibly the corre- 
lation between stature and humerus, while leaving that between stature and femur 
the same. The sensible reduction of correlation between stature and humerus (:809 
to :771), which is found as we pass from male to female, does not arise when we 
select a group of males with their femur and humerus of the same > length, variation, 
and correlation as those of the females. 


ILLUSTRATION ITI.—Jnfluence of a Selection of Stature in Modifying Femur 
and Humerus. 


Let us select a group of French men having the same height and variability in 
height as French women, and calculate the changes which will arise in their femur 
and humerus. 

Here the selection is given by : 


hh = — 1224, ass 5°45. 
We now need only the earlier formule of this memoir. From (1.) we find 


M, = m, + °349h,, M, = my + ‘227 h. 
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These give for our case: 


M, = 40°95, M, = 30-24. 


Thus a group of males, selected to have the same stature as the females, would 
have a slightly shorter femur and a slightly longer humerus. A slightly longer femur 
in woman and a slightly longer humerus in man would thus appear to be sexual 
characters. 


Turning to the variations, these are given by (xlix.). We find: 
o= 236, 5, = 1°53. 


This shows us that while the selection would give the same variability of humerus 
to the men that women have, it would fail to produce the reduction of variability in 
the femur, which is characteristic of the women. 

From (li.) we deduce 


Da = "BUS, Da = "8006, 
while from (lv1.) we have 


Thus we see that very small changes would be made in the correlations, stature 
and femur, stature and humerus, and femur and humerus, if we selected French men 
to have the same size and variability of stature as French women. ‘The explanation 
of this lies in the nearly equal absolute variability of the two sexes with regard to 
stature, for, as we have seen, 14 1s the selection of variability which modifies correlation. 
Looking at the table of values on p. 36, we see that the largest difference of 
variability in the two sexes lies in the femur, and accordingly it is from a selection of 
femur that we should expect the greatest differences in the variability and correlation 
of the two sexes to have arisen, but even this difference alone would not account for 
the observed sexual differences in the correlation. Indeed, 1t would be surprising if it 
did, for the selection of other organs, notably the pelvis, must have played a con- 
siderable part in the differentiation of sex.* 

(9.) I shall now proceed to a series of problems, which will show the application of 
results, such as those obtained in this memoir, to questions which arise in dealing 
with inheritance and selection. If we suppose a general population to have statistical 
“constants,” which remain constant at any rate for a moderate interval, we still want 
to know not only the error which may arise from a random sampling, but also the 
sort of effect which results from our sample being too much drawn from one kind of 
environment, from a rather limited class, or from any other practically necessary or 
unconsciously introduced limitation of the random character of our sample. 


* «Primitive man and woman are more nearly equal in size, variability, and correlation than highly 
civilized man and woman” CG Roy. Soc. Proc., vol. 61, p. 354). 
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ILLUSTRATION [.—To find the Influence on the Intensity of Parental Heredity of 
the Selection of Parents. 


Let the subscripts 1, 2, 3 refer respectively to father, mother, child. Let us first 





select one parent—say, the father-—very stringently, 2.e., s/o; = u is small. Then we 


need only equations (xlix.) and (li.). These give us : 
oe 


ee AS 


and 


The first may be written 


Py. 

E e, E ; ' 

la. = Tap SC? 
er ER EE 1) 
e GN 


1 





or, we see that rı, will decrease, as u, decreases. Thus if 7,, = "A we have for 


2, == 2 132 EE 
py = 1/4, Da = "1085, > 00 = 006: 
fe == 18: Vn = "0049, KE T= ESO, 


It is clear, therefore, that the correlation of parent and child will be much reduced by 
such a selection. On the other hand, the regression coefficient will not be altered, t.e., 
Sol j3/8; = ainsi, as we have seen. Hence in problems of heredity, where we 
suspect a parent to have been highly selected, we should seek for the regression of son 
on parent rather than for the correlation. Thus in the case of Basset Hounds,* some 
if not all the reduction in correlation between sire and offspring may be due to 
selection of the sire. A test of whether the reduction in correlation is due to 
selection of a parent ought to be given by a comparison of a and SZ, We cannot, 
I think, suppose, unless natural selection be very stringent, that o, differs much from 
o, Hence it follows that 3,/s, ought’ generally to be large, if there be selection of 
a parent. We can hardly test this point effectively in the case of the Basset Hounds, 
owing to the nature of the classification. In racehorses, although the sire appears 
to be far more selected than the dam, there is not a great reduction of the coefficient 
of correlation between sire and offspring,f au appears to be less than 3,, but not so 
greatly and certainly less, that we can be surprised that the correlation of sire and 
offspring is not much less than we have found it for material in which selection of the 
father is certainly far less marked. We must accept the warning as to the reduction 


* “Roy. Soc. Proc., vol. 66, p. 157. 
ft See ‘Phil. Trans., A, vol. 195, p. 93. 
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in correlation produced by the stringent selection of one parent, but we must 
remember the complexity of the factors—the variety of other influences at work in 
selecting and modifying selection—before we lay much stress on this source of 
alteration in parental correlation. 

Now let us deal with the case of both parents selected, and suppose their selection 
given by zo = py, S/o, = po, and the change of their coefficient of assortative 
mating from 7, to py We have from formule (lxviii.) and (lxx.) by a little 


rearranging 
Ze = o {1 des Big (1 Ka DM KÉ Bag" (1 DS ul SE (rio ES P 19/41) Bisbo} 


23013 = 0% [913 — (1 = pu) Biz — (ia So Prats) A E (sot) 


where Bis = ("is — agrio) / 1 — 112%), Bas = (Pos — Mis") / (1 — 719). Now let us take 
special cases to bring out points. Let us suppose fri = 0, or no assortative mating to 
exist, and let us enquire what change would be made in parental correlation if we 
selected parents who had assortatively mated, without altering their variability, %.e., 
let us take u = 1, u, = 1, we have at once 


So la EE ee 
Selig = 0 (ris + Pits) e, Uess 
whence we find for ris = Ten = '4 and 

Py = ‘1, 23/03 = 10149, ti = "4399, 

Pig = 2, =e == 2003, 

a= 3, == 150469, = A067, 

ee = 1:0621, = ‘5271, 

Dua = "A. = 10100, = EE 

Pis = I, E oe = ‘6963. 


Thus the general effect of assortative pairing of parents 1s to increase the correlation 
between parent and offspring sensibly, but not to very rapidly increase the variability 
of the offspring. Thus marriages within a class would, if heredity statistics were 
collected for a class, tend to show increased parental correlation. Very high assorta- 
tive mating no doubt occurs with some forms of breeding, and we may well find in 
such cases higher values of the parental heredity than we should obtain for a popula- 
tion of the same species with random mating. I think this may be an effective factor 
in the raising of the parental correlation in the case of coat-colour in thoroughbred 


horses. 
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Now let us see what happens if we select both parents moderately. As test cases, 
let us take q, = u, = H and ‘5, and for the extreme = 0, 
We have at once 

















e “ny 


SS oa 68 eT aoe) a Se EE 
Zei — Fa ¡dp (1 al» Pro) D è H $ D D D (Ixxv1.), 
Henee we deduce 
poy = pry "D p= p= "A, p= pa = 0. | 
ae A 5 e A E 
| 23/03, Tis. | 23/03. tis Z3/03. (GË | 
AA E Se A ls POE Rete 
di -9515 -3700 -8764 -2865 | -8246 0 
“2 -9622 -3991 -8809 -3093 | -8246 0 
3 "9727 Ev ‘8854 ‘3316 | "8246 0 | 
d :9832 4556 8899 3535 | 8246 0 i 
5 -9936 -4831 ‘8944 -3750 8246 0 
A Ei | 10438 -6131 9165 ‘4762 | +846 o | 
| | | | 


This table is very instructive. It shows us that selection and assortative mating 
are factors of opposite influence; that if selection be only moderate, then with 
considerable assortative mating the coefficient of parental correlation may be 
increased, but if selection be stringent, then assortative mating cannot counteract, 
even if as high as "D. its destructive influence on parental correlation. 

For example, if we take parents remarkable for some intellectual or physical 
character, say with a variability only a very small fraction of that of the general 
population, then, however proportionately we might pair them, we should find their 
relationship to their children, as measured by the coefficient of correlation, very 
sensibly reduced below that of the general population. 1 think we have here the 
reason why Mr. GALToN's Family Data, which were drawn from a rather narrow 
class, and had only a small coefficient of assortative mating, give so much smaller 
parental correlation than my own Family Data, which seem to me drawn from a 
wider class, and have a considerably higher assortative mating.* 

It will be clear that with factors like assortative mating, natural selection, artificial 
selection of breeders, unconscious selection of material from one class or one environ- 
ment, modifying our coefficients of heredity in one or another direction, we can hardly 
hope for more in practical statistics than an approximation to the strength of the 


pure inheritance factor by dealing with the average of as many races and characters 
as possible. 


* The work for Mr, GAvron’s Family Data is given, ‘Phil. Trans., A, vol. 187, p. 270. My own 
results are as yet unpublished. The average value is about +45, as compared with Mr, Garron’s +34 
VOL. COA. G 
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- ILLUSTRATION 11.-—To Find the Influence of Parental Selection on Modifying 
Fraternal Correlation. 


Let the subscripts 1, 2 represent the parents, and 3 and 4 two of their offspring. 
Let us first select one parent only, the selection being given as before by side: = py. 
Our formule will now be (liv.) and (lvi.). So far as the change in variability is 
concerned, we have already discussed it under our first illustration, so we need only 


consider : 
r. (Poy — riga) QU SS Di + py” "Peg SET, 
GR J — A mre il — 0 — ar) 


Now ris = rip if the offspring are of one sex; hence: 

















(zen) 


Meal — paí) ue 
EE ee 


1 Sr SOL — pa?) 
If we take 7,3 = '4 and fu = °5 as reasonable values, we have 


Er 16p, 


3t 84 1 ën (lxxix.). 


Ta 

Thus r}, will be greatest when p, 1s greatest, 2.€., when there is no selection, and 
will decrease with pı until it reaches ‘4048, when p, = 0, or there is selection of 
fathers of one value of the character only.* 

The selection of one parent only does not, therefore, immensely modify the 
correlation of brothers. Still, if we work sensibly with one class of the community— 
say, men of genius—we should expect to find their sons rather less like each other 
than if we worked with the general population of brothers. 

Now let us select both parents. Here again the variability of the offspring has 
already been dealt with. We are concerned with equation (Ixix.), and we shall put 
Tyg = Ty = Tog = To, = Y, Or make parental influence equipotent for the two sexes. 


Hence 


Biz = Baz = = Br, = Ss Boy = p = => 


where 7 is the parental correlation, and 7,, the coefficient of assortative mating. 


Hence we find 


L+ e 








Pg BIL — m? +1 — us + 2 Une — Prot) $ doe A). 


2 
Toy == o 
í 1— — py,’ + 1— pe + 2 (1% — PrigpMa) $ 


u 
To reduce to numbers, suppose q, = po, and 7. = 0 for the general population. We 
have 


_ Ty — 2r? {1 ar (1 + pig)t l foe te Ga) 


To, — 
1 = 2 11— 2 (1 + py} 








D 
KR e 
* Tn general the value of rs, ranges from sg, down to SC S T 9 
— 113 
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Hence if we put 75, = A, and r = ‘4, 


e —= 34 32%" (1 pel 
34" 68 + 324 (L + pio) 





(lxxxii.). 


The following table will suffice to indicate the changes which take place, when we 
give a series of values to u, and ps, Thus the first row gives the influence of selecting 
parents without any assortative mating. We see that with increasing stringency of 
selection the reduction of correlation is very considerable, and that with such selection 
the influence of assortative mating becomes less and less. Nevertheless, assortative 
mating can produce quite sensible results, if there be little or no selection. I am, 
indeed, inclined to think that a good deal of the high values found for the fraternal 
colour correlation in the thoroughbred foals* is due to much assortative colour mating 
in sire and dam. Of course 1t cannot be all due to this source. 


VALUES of Fraternal Correlation with Parental Selection. 




















| 
eae meer e EE p="4, E AR =O: 
E + o | OE T E Ba SERENA 
| pas -5000 "4349 | la “o 169 -2783 -2647 
(p=? ‘D155 44TT — | + 3802 -3209 -2796 -2647 
mass "2 -5301 4599 : 3889 -3256 :2809 2647 
NEE ‘5438 -4716 -3974 -3303 -2823 2647 
mass H -5690 "4935 4137 - 3392 284.9 | 2647 
o pia=1 «6912 | 7514 | 4508 - 3609 “2914 | 2647 
| | V 




















On the whole, I think, we may conclude, so far as the relative influences of sexual 
selection in the form of assortative mating and natural selection go, that : 

Both sexual and natural selection can sensibly modify the intensity of inheritance 
as measured by the coefficient of correlation, the former tends to raise, the latter to 
lower, its intensity. But the effect of the latter, if at all stringent, ts to completely 
mask the effect of the former. 


Tn tact, we may write 
ee 


ro = ]— : 
: 1 = 2? + 2%? (L + py) 


Hence the smaller p,°(1 + pip) the smaller will be fraternal correlation. This 
varies as the square of m and only as the linear power of 1+ pı} Thus we see at once 
why stringency of selection 1s far more potent than assortative mating. | 

ILLUSTRATION IlI.—To find the influence of selecting two organs A and B in 
a parent, on the correlation of the like organs A’ and B' in the offspring. 

Let the organs in the parent be denoted by | and 2, and in the offspring by 3 and 4. 
Suppose the organic correlation of the two organs in the general population to be +”, 


* “Phil. Trans., A, vol. 195, p. 93. 
G 2 
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so that r = rs, = 1” before any selection takes place. Let r, the correlation of the 
organs in the parent and offspring be supposed to be the same for both organs ; then 
r= fi = 7, Finally we have the coefficients of cross-heredity, zu and 73. These 
must vanish if there be no heredity and no organic correlation, and should be perfect 
if both these are perfect. Hence we will take 1,, = r, = 1” as a probable hypothesis.* 
With these values of the correlation coeficients we easily find 


Du = r = Boy Bag = By = 0. 
Hence from (sen, 1 and ([xix.) we have : 
2 =o E re, SS eh ew ers)... aged), 


eee PLS) + pe eee 
A a ee h 











(Ixxxtv.). 


For simplicity, suppose the stringency of the selection to be the same for both 
organs, then: 
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Acel 
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E oe 
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yee X Een (Ixxxv.). 
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The following table indicates the value of t/r" —- 








| | | 
| Ui zs lL. UE e SE | m pnd. | 0 
EE o e a Md o la E E E A A ee 
| | l ' | l | | 
y =1 | i 1 | I | i | 1 1 
y=0:8 | -9680 | oan | ‘9872 |  *9941 9985 1 
y=0:6 | +9860 9565 | “9743 | on | -9970 | 
y=0:4 | -9040 -9348-9615 | -9823 | -9955 1 
y=0:2 8720 | 9181 | “9487 | +9763 9940 1 
y=0°0 | -8400 8913 (mom | moi | 9924 1 





It will be clear from this table that if the selection be at all stringent, no reduction 
of organic correlation in the parents will affect substantially the organic correlation in 
the offspring. 

On the other hand, if y be > 1, we can have considerable modifications in the 
value of the correlation, even if the selection be stringent. 


* See ‘Roy. Soc. Proc., vol, 62, p. 411. Ihave a good deal of data on the value of these cross-heredity 
correlations now reduced and soon to be published. 
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Thus we have the following values of 1,,/7’, if :— 





Ae aii ine e EE a EE A e a a A s wé AAA O RO E ENER EEN EEN ER ENEE at S ati E mt bt rn 








pl: p= 8. EMA payaso 4. EE EES 
EE EE 
y= 2 | 1:1600 1:1087 1:0642 1:0296 1:0076 1 
= 10 2: 4.4.00 1:9779 1:5775 1:2662 1-0681 1 
= 50 8:8400 6:3243 4:1448 2:4492 1:3705 1 





Lastly, if ui = 0: 
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Va Pi E 
Ee 
mens 4 fy” 
OF, ts;/Pi 18 given by : 
EE EE LEE 
| menk Mee E pee 0. ere | jg. py = 0 
ERENNERT EEGEN EECH ke 
Y31/p12 "1600 | *1087 "0642 "0296 "0076 d 
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Thus, even 1f there were no correlation between the organs A and B in the general 
population, still a selection of parents in which such organs were correlated would 
lead to offspring with correlated organs A’ and D The amount of such correlation 
would only be °1600, if the variability of the parent were not selected, and would 
diminish rapidly with stringent selection of variability. Still :1600 is quite sensible, 
and would, if the selection continued for a few generations, continue to increase. 
Thus we see how selection of a pair of organs in a parent may increase or even create 
correlation between the like organs in the offspring. 

The reader will find other interesting illustrations in tracing the influence of an 
absolute selection of one parent only on the correlation of the offspring, e.g., relation 
between pairs of foals which all have a common sire, the influence of selecting an organ 
A in the sire and an organ B in the dam on the correlation of the organs A and B in 
the offspring, the influence on assortative mating of selecting parents of men of 
genius,” and in many other problems. 

(10.) It is not without value to consider how pj, arises in the case of natural or 
artificial selection. Suppose we have two organs, A and B, then we shall endeavour 


(i.) to give these definite values, say a and y, but we shall not be able to get all our 
individuals with such absolute values, we shall select with certain deviations from 
wand y, given by w=a+a' and y=y-+y’, say. Further, we shall endeavour to 

* This is a peculiarly interesting case, for if we select men of remarkable intellectual ability, we should 


expect to find both parents above the average of the general population, but with a negative correlation 
between them amounting at a maximum to — ‘1905, 
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: De e SE / = ae pp ee 
make y some function of x, yy Ef ol ory +y = f(a) + x f a) + E J” (ae) +... 
by TayLor's theorem. But y =f (x) and al is small, so that our attempted relation 
will be of the form : 


} F 
E 


Here m is the slope of a line to which we endeavour to confine the selected organs. 
But we shall not be quite able to exactly hit this relation off; when a’ = e, we shall 
find that y = me + y, where y is a small unavoidable error in selection of A uncor- 
related with e Thus, if s, and s, be the selected variabilities, we shall have : 


TET 1 ; 1 9 
Ss Oa e : S(me+ yy, ai = S S (€). 


Oé 


1 
Therefore : Ss == Sm). 
2 EE EE 
S (yx Sí(me + me ms? ms 
E Ph = Syw) _ Sim DEl EE, 
7 NS\So NSIS» Dies SS 
d 
and therefore : =D 53. DE 


Or, pa is at once obtained from the slope of the line m, by which we endeavour to 
fix the relationship of the organs A and B. Or, again, we may look upon 
SE — pi? as a quantity measuring the mean divergence of the B organ from that 
absolute fulfilment of the relationship between A and B which we are striving to 
attain. Thus p is a quantity which naturally arises in every attempt, whether 
artificial or natural, to select organs having a definite relationship to each other. 

Much the same considerations arise when we select three or more organs. In each 
case the selected coefficients of correlation are constants which enable us to express 
(i.) to a first approximation the form of relationship we are aiming at, and (1.) the 
average degree of divergence from absolute fulfilment of this relationship. 

Thus, without regard to any particular distribution of frequency, the s’s and the 
ps are the appropriate constants to express approximately the nature of any form of 
natural or artificial selection. 


(11.) On the Probability of Survival and the Surface of Survival Rates. 


In the course of the present paper I have assumed that when measurements are 
made on any population for a complex of n organs, the frequency surface may be 
taken as approximately normal. If this holds for the population before and after 
selection, and measurements made on many groups at different periods of life seem 
to indicate that it is approximately true, it follows that we can determine the form 
of the probability of survival as a function of the means, variations, and correlations 
of the selected and unselected populations. 
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Let the unselected population be given by 
Z = Zo expt. — Semi + cyto? A... F Cain 
+ Geste HF... + Bu sl Gul, o (Uess) 
Let the probability of survival be given by 
pP = Paf (&ı— ki, ko Ly — Ks, o. By kh) . . (lxxxviL), 
where fis at present an unknown function, which is to be a maximum for 
El Et, a he ie. 
and, if the selection be at all stringent, to take rapidly decreasing values as 
Ly — ky Xy — ko Xa — ko, o o. La — k, 


take increasing large negative or positive values. It will be clear then that the 
individuals who are “ fittest to survive,” ie, have the smallest death-rate, are those 
whose organs are defined by : 


dr _— k,, Lo — e ss o Ly = A 


and fitness generally will be measured by the closeness of the individual to these 
“ fittest” individuals. 

In order to find the surface of survivors, immediately after the selection if growth 
be taking place,* or at any later stage 1f growth have ceased, we have only to multiply 
Z by p, or: | 

EE be ee en ee lee 


is what in the earlier part of this memoir 1 have termed the selection surface. Now 
if this selection surface be itself normal, it will be of the form: 


z = ën expt. — $ {bn (2, — hy? + dog (2 — ho +... 
Se Din (Ly a li)? SC 2019 Los GE b,) (ke SC bo) 
Or (Ea — Ani) (£a — hn) |. (Ixxxix.). 
Here, as in the value of Z, all the constants b,,, by, . > Dan: Dia... Onin are known 
in terms of the variations and correlations. If there be selection of q organs only 
out of the n, then Dana, + Dum Olgi Doan, + + Oaia Will all be zero. Since by 


Equation (Ixxxvill.) p = 2/Z, it follows that the function f which defines the prob- 
ability of survival must be of the normal exponential type, or 


* T propose to deal in another memoir with the important problems of slow selection during rapid 
growth, and of secular selection during several generations, 
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fla, — Eu, Zy m Aa, Lym ul = expt. — 4 1a (0, — EP + gg (Hy — By)? +. 
Y Sos (20, — Ey Lg — By) E Llao (na e hna) (a — Al e Ee). 


Thus, to determine the probability of survival, we require to know the values of 
the oa and kam terms of the Ga Xs, and en The shortest method of finding p, Is 
to put ms hy, a = Ba, % = Bn ahd then note that : 


ES € F dy 
El, EE Än, oe Uy = = A 
TE ee XP SRIOID EEN ` ENEE te mee a h aware ot ge 


E SE SS wn O 
£0 Vd (2, = hy, da == kt, coe hy == ih) ( ) 
Since pegs, ml 2 ae", 


we see that the relations for p, given z and Z, and for z, given p and Z, are cyclicly 
interchangeable if at the same time we change on to — Cip Cag t0 — Cog... Cy to 
e Cige e Cua TO — Cuco Lf o, 07... 0, be the standard deviations of the 
unselected population, this amounts to changing o, op... 0, to Y— lo, /— ley, 

., — lo, respectively. Thus the results which give the probability of survival 
in terms of unselected and selected populations can always by an easy interchange 
be used to obtain the selected population from a knowledge of the unselected 
population and of the probability of survival. 

Let the unselected population be defined by M,, Mas. o. Mu Oy Ox + + + Ty, and 72, 
Tiss Togs o o e Pnl m | 

Let the selected population be defined by m, + Hy, mg + Ha, ... ma + Ha, zz 

E E E EE E 

Let the constants of the probability of survival function, or Ou, Gy... © Qan 
Cía + + » Gaza be expressed by Sp Sa, + + So Plas Pig + > + Pronn AS if ib were a normal 
correlation surface.” 


Then the problem will be solved, if we know the a, sn, and oa in terms of the ds, 
va, H’s, ¥s, and ra 


E uatin the sc uares products and linear terms in the XS In the e uation 
H I 4 H 
p= z/Z, we ha ve at once the Sy stem ` 


e 3 
yy Gs Den == Cons i fad on ate 
ee ee Kee 
(un Some CR wee Ca: | 


for all values of u and v from 1 to n. 


E ll, oe Codes A eg nal. A eege Cry ke, 
SSES ba H, 153 GH. KEE eg bah, E ie? EE 


¿ Cin), 


for all values of v from 1 to 7, 


ate 


* These must not be confused with the sj, $», .. . Sy Pi» Piz P23) ++» Pag- Constants of the y selected 
organs of the previous discussion, The new quantities may be in part imaginary. 


SELECTION ON THE VARIABILITY AND CORRELATION OF ORGANS. AY 


If we now substitute from the first two equations for the og we find : 


(6, on D ky + (En E D ky + coe SR Dës Lee D hy + e.: + (e SR Din) ky, 
= — baH, — beH, —...—-b,H, ...—0,.H,. . (xciv.). 


Now let A be the determinant 


Cp Dun, Ee Dir Ce = Dis, hae ven Cee Dy, 


e)? 


Ca — Dar, Gem Daz, Cog. E Das, E . Cay — Don 
gt Da, C53 — Dan Car Dan, Can oo 0 
| 


Uy Dar, Cu? D C3 TRR ‘ . . . Cra D 


ged 


and A,, the minor of its uv™ constituent. 
Then we have at once: 


A X ko, = Any (H, + b,.H, + s. + D.H. 
+ Av (baH + byH, +... + b2,H,) 
+ As (ba H + bel, + basHg. o + DH. 


SS Aen (ba H, + (SE, + o. + ele Ee 


A X k, = Hy (ban + bit + digA Lk... Oy Ann) 
Ss H, (0314, SR Gah, Sk DogAns eee Boyden) 
+ H, (bz1An + Base + DegArs +... + bah) 
+ H, (OnAn + due + du Ar... BOA) . . . (ev 


Thus k, is determined in terms of H, H,,... H, which define the maximum 
frequency of survival. 
In a similar manner by making the proper changes indicated above we have: 


A’ x H, = ki (ay Al + As bys Aig +... + An) 
+ k, Lan A + dre F AA ho... + Gah aal 
+ k (UnA n F aA wn F AA Lt... E Any den) 
ey (anA e aA E ag ig me Ei in) oe (ev), 


where A! is the determinant 
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(Cz Se Goal, (Gy, ae Cts), (Cs Se Asa) o. ( 
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(Gat Aa), (C + arh (Cn Lal... Iron F am) 
and A’,, the minor of its uo constituent. 

Now it will be clear from these results that as a general rule it is impossible for E, to 
be equal to H. In other words: The individual most frequently met with in any given 
selected community, 1.e., the mediocre indimidual, ts not the individual fittest to survive. 

It is only in the limiting case of natural selection being so stringent that one type 
of individual alone is able to survive, that the fittest class has a numerical majority 
over any other class of the community. This seems to me an important, algebraically 
almost self-obvious truth, and yet one which is very much obscured by the use of such 
a phrase as the “survival of the fittest.” 

Of course, if there be continuous selection, or an environment so stable that the 
probability of survival remains constant for a long period, there will be a gradual 
approach, never theoretically an actual identification of the mediocre and the fittest. 
But in actual nature the environment, at any rate so far as it depends on climato- 
logical conditions, must have a long period as compared with the vital and reproductive 
periods of innumerable forms of life. A hard winter, a drought, a flood, a famine, 
a plague or epidemic of any kind, even if fairly stringent, will rarely, if ever, render 
the most frequently surviving individual identical with the individual who is fittest 
to survive.* Still less will this identity take place in the many processes of artificial 
selection, which are becoming and will more and more become valuable laboratory 
aids in our appreciation of the action of natural selection. The divergence between 
the most frequently surviving and the fittest individual is measured by the above 
formulee for the /'s in terms of the Har 

To complete the solution, the oa must be found from the equations of type 


Ou = Dom — Co, and then from the a’s the s’s and ps follow by the well-known 
determinants for multiple correlation : see our Equations (x1.) and (x11.). 

Throughout the earlier part of this memoir I have used only the surface of selection, 
but the above investigation will enable us whenever desired to replace 1t by the 
probability of survival. L will illustrate this by obtaining the formulee suitable to the 


simpler cases. 
* We badly want a name for the selection which acts for a short time and rapidly modifies the adult 
population. It is practically the type of selection considered in this paper. It is epidemic or catastrophic 


in character. 
+ The point is of considerable importance, for more than one influential writer has spoken of the result 


of natural selection as the preservation of the type the mortality of which is /east under the given 


conditions. 
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(12.) Case (1.).—Selection of a Single Organ only. 
The original population is given by 


Na — Lie FE, 
and the curve of survivors by 
z = z e` lr HEES 
0 . 
The probability of survival 1s | 
p — Po E Ge 
where we easily find, if 3/0 =A, 


SEV Ga. Ga ON 
s= Se A CERN (xevill.), 
E eg 
and | a= AE Ee, cea): 
0 


As an illustration consider a selection from modern French peasants, which should 
reduce the mean and variability of their cephalic index to those of the Libyan race. 
French peasants :— 
mos 19786, o = 3841. 


Libyans :— 

m + H = 72938, Dames iy 
Hence : H = — 6'848, Ai ol 
These give : CSS dale, LEESCH 


Thus for such a change as 7 points in the cephalic index to take place by selection* 
we should have to make the “fittest to survive” of such a ridiculously low cephalic 
index as 64'074, and such a high variation as 4'370. 7 


We find Po =-51:0474 n/N, 
and accordingly the probability of survival given by 
p — 51:0474 = o — (e+ 15°712)21(38-1939), 


where N are the number of Frenchmen converted into n Libyans so far as cephalic 
index 1s concerned. 

I have purposely taken a somewhat extreme case of selection in order to illustrate 
how widely the most frequently surviving individual can diverge from the fittest. 

In this case, if the chances of survival (1.) of the fittest, (11.) of the individuals most 
frequent after selection, and (111.) of the individuals most frequent before selection, be 
C,, Co, and C, respectively, we have : 


* This is, of course, supposing the change to occur by catastrophic selection and not by a continuous 
secular selection, see footnote preceding page. 


H 2 
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or, the chances of survival of an individual of the fittest type would be about eight times 
as great as those of an individual of the most frequent type after selection and about 
700 times as great as those of an individual of the most frequent type before 
selection. If »,, Ge vz be the numbers after selection in the three classes* of the 
fittest to survive, the most frequent after selection and the most frequent before 
selection, we find 

Peay Pe 00802 vel “00076: 


In other words, the most numerous type before selection 1s still after selection about 
6'7 times as numerous as the type with the least mortality, and this latter type is 
only about yr as numerous as the type to be most frequently met with after 
selection has taken place. 

Thus, although there would have been a very great evolution in cephalic index, 
due to a fairly stringent selection, the fittest to survive would always have formed 
but a small fraction of the dominant type. Even if we were to replace the selection 
here considered by a gradual evolution spread over several generations, we should 
still reach in the main the same conclusion, 2.e., that natural selection never proceeds 
by the survival of the fittest, or the survival of those with the least death-rate. These 
will always remain a small fraction of the community—they are the goal, but often 
the very distant goal, to which selection tends to shift the population. 


(13.) Case (Di. )—Selection of Two Organs. 
In this case let the surface of survivors be: 


Ge a (a, — H. Bhi (e a —- H,) (2 — H,) — H,) (2, — H." k 
>, 20) expt. — 9 E SE ta 2) >a SS 1,2) SR 5? (1 Ke t) o e) 


the original population : 











S ao oy EE o H S 
Z = Ly expt. — $ T a a er a o. (cL), 


oy" (1 — 73") 010 (Loe ap) Gh a). 
and the curve of probability of survival : 


Ro — Au Zë Loi — Ay) (a — Ao) Ug am Ba)? hy)" — db 
) = Py expt, — + bl oS eee ee ae GC 
P = Po ent el aL == Pig 2) 8,5 (1 — Py) Se “(1 — Pig”) ) 








Since : KEE 


* By individuals of a type or class is meant here, as elsewhere in this section, all the group falling 
within some small definite range of variation lying round a particular value of the organ (eg, m, m + H, 
or m +), which defines the type or class, 
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we find at once 





























1 1 1 Ee 
GE AE e e "EE 
BE = pi) En ET AL) í 
1 1 1 1 SH 
EEN E EE EE ee E 

SEL Se E EE a E \ 
| Pia — E Lu Bé CAE (ov \ 
53, 1 — Pr 227 (1 — 2) oalr) ` ffe E 
Ges Ay es Daf ` BEE EE Hato de (ovi ) 
E SE 5,5, (1 — fra”) On eS oe a) SiS O ES 
Pr LA E 7 a ea Y "fl, ` H, (evi? ) 
35 (1 — py) "32 (1 — A) Zästl të * 32 (1 — ns) ES 


Let > Ve = As 2/0, =A, measure the stringency of the selection, and 


— Tr a 
= /1= 7 E measure the change 1n correlation, * Then solving the above equations 
me 
we find : 
Ya Tat 22 ¡A 








SE eg , Get, ), 
Ca" Y lr? SO i.) 
EE (cix.) 
A + (1 — tia?) A 227? + Zare, Wi D 
B OROEN EE 
Sa oe ( E S (Gia); 


1 "ail ig a A tL he TAr Aa H PUAA 


ley le 7% — Ag? + hM 4 


5) = 0, D E 


NA GN EE i 
e Sg l l i f : (exi), 


Ko — ly hAg — "Ai T h, 1 = Ty Af A hahh (exii ) 
G3 2 B e e BS E 


Boa Lr yg? = A Ag EU 01) APA + 241991242, 








where 
Similarly, if the original population and the curve of probability of surviving 
or of survival rates be given, we have to find the selected population : 


Pro + Pig” Ky Ky 


Lt Pee Sl + PK, 





e (exiii.), 


* Tf ro = cos d, tig = cos D, p = sin D/sind. The quantity D has been conveniently termed the 
“divergence” by Mr. SHEPPARD. Hence p is the ratio of the sines of the selected and unselected 
divergences. The above formula for pj2 can be at once changed into one suitable for trigonometrical 
logarithmic calculation. Let sin a = pàn sin d = pàs and pra = cos ô; then, if A be the side of the 
spherical triangle, of which «1, a» are the other sides and 6 the included angle : 

sin 43 = Ne de 5(D — A) sin $(D +A) 


COS 21 COS La 
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har L= p a and eae 
where ye ioe Er S/o, and ky = Salz 


A ct A < c<cF——— A AA (cx1v.), 
o H — ` 2, eye 








Ns MENA + Cm By?) E: 


2 ee Sa a NO ee (cxv.), 
(i Sig be oP hy a (1 — fig") "Kg" — 2h iQ? Qk Ky 





o o Y 07 Y 


hi, TK — Pik E 
E 1 puta : 
ere! TS i er (cxv1.), 


3 es js e DH 3 9 A 4 
H, bh Tia — Pigia ep hy L — Pig + ef — Pigia 


O 9 On y o D y 


be 








(exvii.), 


(14.) Fllustration.—The following results are taken from the paper by Miss ALicu 
Ler and myself already cited :—- 


French g. Aino d. 
Femur . . . . m, = 45'228 centims. 40'770 centims. = m, + H. 
e cs 1:898 e O 
Humerus . . .m,=33010 , 29502 | cs = Mo, + H, 
os = 1538 , 13483, = ZA 
Correlation of? i | 
Femur and e, . igs “8421 "8584 = Uy 


| 
CED | 
Humerus J) 


As indicated by the symbols above, let us select from the French a population 
having the same femur and humerus relations as the Aino. 
We have at once: 
H, = — 4'458, H, = — 3508, 
e 80017, M= 87321, 
u? = 9047, 


Il 


Whence we find : 


| 


Pr = ‘9027, 5, 3:4870, 5, = 28786, 
kı = — 178464, ky = — 158547, 


Tb 
Po = 208,425 y; » 


where n are the number of Ainos which can be obtained from N Frenchmen. 
We have accordingly the following form for the surface of fitness to survive : 
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qo 
wor 


had 


p = 208,425 À expt. — 4 {222056 (æ -+ 17:8464)? + ‘326,960 (y + 15'8547)? 
| — "486,450 GL 17:8464) (y + 15°8547)}. 


Now it is clear that if we wanted by a “catastrophic” selection to convert the 
French into something resembling the Aino, we should have to give the least death- 
rate to those French with femur corresponding to m, + k, and humerus to m, + kə, 
or to the dwarfs with femur = 27:382 centims. and humerus = 17:155 centims. ! 
By no other means could we shift the modal value of the French population down as 
low as the Aino modal value. The physical meaning of this is that we have been 
compelled to put on an excessive death-rate for the bigger Frenchmen. 

An interesting point of our work 1s that 


k = 12514 H, + 3:4971 H. 
k, = — 1:0266 H, + 58242 H. 


whence we see that while a selective reduction of humerus is far more effective in 
reducing both femur and humerus centres of survival than a reduction of femur, a 
selective reduction of femur occurring contemporaneously with that of the humerus 
actually tends to raise the centre of the humerus, ze, the coefficient of H, is 
negatwe. 

Now let us consider the frequency of survivors per unit length, say centimetre of 
femur and humerus, at different points. The surface of survivors, je, the Aino 
population, 1s 


eng 
2 





(nit 3 2120 Zei 


1 {GP 2r Ge - Hp) (y - Ho) (y — Hey) 
Ki) See RER e el $ 
be? 


a. TT II TT e 
2172129 es thy”) 


If we put x = 0, y = 0 we have the frequency after selection of the original 
population type ; if we put « = H,, y = H, we have the frequency after selection of 
the new population type ; and if we put « = ki, y = ky, we shall have the frequency 
after selection of those best fitted to survive. If these frequencies be va va, vy 
respectively, we find on substituting the numerical values that 


Dios Van tae Qe? le 082 280. 


Thus the most frequent type of the new population is now about thirty times as 
frequent as the old most frequent type, while the type most fitted to survive has 
practically no existence at all. It probably lies outside the actual boundary of the 
French population. 

Here really arises the question as to how we are, in any actual problem, to fix the 
ratio of n to N, or, what amounts to the same thing, to fix a practical boundary to a 
given population. Such a boundary must be conventional, but I think that for 
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practical purposes we are quite safe if we assume that an individual who occurs only 
once per thousand can produce no effect on the physical evolution of the population 
as a whole. 

Now the form of a correlation-surface for two organs, x and y, is 


N SC 1 gar tray Y? 
e 1 op? to? TiTa od 





e : 1 ae Zomm 2 ; 
Let us write 2 = 7 s ” E GE then x = a constant gives a series of 
1— lo? oo, | off’ 2 


similar ellipses which are the contour lines of the surface, or lines of equal frequency, 
ie., giving individuals with equal probability of occurrence. Let the equation to 
these contour lines referred to their principal axes be 


a X? y 
E ee 
Then we have at once : 
i 1 = 1 P 1 a? + 07 
Ai TD (1 oeël--aä1 1? 00, " 
1 1 qa 1 





AB T opos (1 -- 1°) K TGs (ee Py = (ies ECH 20," 
or, AB =o vl — r’, A? + B [= o HH o. 


Further, if db be the angle the A principal axis makes with the axis of x, we 
have :* 


These fully determine the principal axes of the frequency surface. Now consider 
the frequency between the elliptic cylinders corresponding to « and x + òk; we 
have tt 

= z X 2rABkdr = 2 X 2T V 1 — 1° nde = Ne kdr. 


Hence, if N, be the frequency outside any contour k, 
a y Ge 1.9 < Pre List (ET? 
N sehN ois sc: Nei, . . . . . . (xviiL). 
e K 
For N, to be gg of N we have simply 
KS > whence x = 3:716,923. 


* For easy calculation put y = Joy? + 0%, tan y = 09/01. Then we have at once if r = cos D: 
A = y COs y, B= y sin x, 
where : sin 2x = sin 24 sin D, tan 2h = tan 24 cos D. 
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This will enable us to determine our conventional boundary to effective population. 

Now let us refer our non-selected and selected populations to their centres and 
principal axes. ` 

We find for the contour curves :— 








Unselected Population | Selected Population Surface of Survival 
(Ge, French). (Ge, Aino). (Ge, Rate of Survival). 

x centre (femur)... . 45 +228 40770 27 :382 
y centre (humerus) . . . 83'010 29:502 17:155 
tan ø (slope tov). . . "601,3775 "670,1454 "807,3371 

PF ae 2:1338 2°2514 4:4115 
Principal axes { p TE 7197 - 5808 9775 
1 in 1000 Dk? SE 10:1615 8+ 3682 -— 


al KR 2:6750 2: 1588 — 


A A A rn enfin ee a m e e een EEN EE EE EES EE ENEE eer 


Referred to its principal axes, the rate of survival is now 


x? ye 4 


ga 
E E 
P= 208,425, DEn en 


Suppose we require to get at least 1000 Aino out of the French population, N, then 
n = 1000. Now suppose the Aino limiting ellipse drawn, then the French population 
must be sufficiently large to give the individuals inside this ellipse. Now p gets 
smaller as we go further from the centre of the survival surface. Hence the contour 
line of the survival surface corresponding to p = 1 must be touched externally by the 
limiting contour of the Aino population, in order that we may get at least 1000 Aino 
out of N Frenchmen. Now, by a graphical construction, I find the major axis of the 
elliptic contour line of the survival surface which touches the Aino limiting ellipse, is 
about 11:44, This gives for the parameter «, of this ellipse, k) X 4:1149 = 11°44, or 
vu = 2°5932. Whence: 


Te: GE 
Pig de 208,425 SNO eTe 5932)" 


gives the greatest possible value of p and the least possible value of N. Numerically 
this gives us N = 7,200,000 about, or we should want more than 7,000,000 of 
Frenchmen to obtain our 1000 Aino by a catastrophic selection. The actual bounding 
contour line of this least possible number of Frenchmen* has for its major axis 
15'285 centims., and it touches the Aino limiting ellipse at the point where it is 
touched by the survival contour p = 1. 

Now let us turn the problem round and ask what is the least population of Aino 


* The least possible to reproduce the Aino, as far as femur and humerus are concerned, to 1 in a 1000 
of the population. 


VOL: COM T 


90 PROFESSOR K. PEARSON ON THE INFLUENCE OF NATURAL 


from which we could produce 1000 Frenchmen by a catastrophic selection. In this 
case the surface of survivals is simply obtained by inverting p, and, if N’ be the 
number of Aino, n = 1000 = required number of Frenchmen, we have : 

À 1 OM 


x y? 
Se al da” SE 


P = 208.425 N 


e 


Here p’ gets larger as we go away from the centre of the surface of survivals, 
and we must therefore make the French limiting ellipse just touched internally by 
the contour line of the surface of survivors for which p’ = 1. The major axis of this 
contour line for p’ = 1 was found by a graphical process to be about 33:91. This 
gives for K, 


Kal (44S) = 8891, a a 10807. 


1000 1 o2 76867) 
3 


208,425 N’ 





7 
Thus: Dl 


leads to N’ = 32,460,000,000 about, or we should want upwards of 32,460,000,000 of 
Aino to produce the 1000 Frenchmen. The bounding contour line of this number of 
Aino has a major axis of 15'890 centims. about, and touches the French limiting 
contour in the point in which it is touched by the p’ = 1 contour of the surface of 
survivals. 

Now the difference between these two unselected populations is very great. We 
see that to get the Aino a very great number of Frenchmen would have to be 
exterminated, about 7000 for each Aino selected; but to get the Frenchmen from 
the Aino an appalling number of Aino would have to be destroyed, upwards of 
32,000,000 for each Frenchman selected. Even if the selection were not catastrophic 
but spread over centuries and centuries, we must recognise what a large consumption 
of life there must be—andividuals destroyed without progeny*—if we are to suppose 
any highly civilized race like the French produced by selection from an apparently 
primitive type like the Aino. Indeed, the return journey in this case seems much 
easier than the upward ascent. Beyond all this we have only made French and 
Aino alike for two organs, and only for one character of each of them! Allowing for 
our conventional limit to the population, allowing for the fact that our Aino data are 
drawn from a very limited population of remarkably small variability, 14 seems very 
improbable that the French have ever been produced by selection from a primitive 
race at all resembling the Aino. The fact that the Aino could be so much more 
easily obtained by selection from the French seems to indicate that they are rather 


* Of course, with a secular selection spread over many generations, it is largely the potentiality and not 
the actuality of life which is destroyed. Still, while the gross number killed among a small primitive 
community may not be large, the death-rate must still be immense. I hope to return to these points 
when dealing with secular selection as distinguished from catastrophic selection. 
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some degenerate offshoot of a race superior to themselves than a sample of. the 
primitive people from which the Circassian races may be supposed to have sprung. 

The whole of this discussion 1s, of course, very hypothetical; no stress whatever 1s 
to be laid upon it except as an Illustration of method, and a rough appreciation of the 
vast amount of elimination which must be necessary to evolve one race from a second 
in the case of organs which we know by measurement to have continuity of variation, 
and: only saltatory changes in pathological cases, which have, as far as we can judge, 
no influence on the mass-evolution which has produced the local races of man. 

But given fair samples of material our method will enable us to determine whether 
a race A——for of course a limited number of characters—could with less destruction be 
deduced from a race B, than the race B from A. It will not therefore follow that the 
path of least selection is that which necessarily was used by Nature. Possibly both 
A and B have been reached by far less expenditure of material from C. Still it is 
something definite in the midst of our gropings after truth in problems of descent to 
have even a rough appreciation of the amount of selective destruction which would 
arise from alternative suggestions. That is why this special numerical illustration of 
the surface of survival has been given. 

The reader will possibly find the matter rendered somewhat clearer by the diagram. 
The femur .is measured along the horizontal and the humerus along the vertical. 
A is the type or mean femur-humerus of the Aino population. Within in the 
continuous ellipse round A the whole Aino population up to 1 in 1000 would fall. 
F is the type of the French population and the continuous ellipse round F gives the 
area within which up to 1 in 1000 of the French population fall. Since the diagram is 
drawn to centimetres of the bones, it will be seen how very small are the limits of 
variation within both populations. P is the centre of the surface of survivals; for the 
selection of Aino from French it makes the “ fittest to survive.” In the case of the 
selection of the French from the Aino, P is no longer the centre of fitness, but the 
“centre of unfitness”; the Aino are killed off with an intensity which increases the 
closer we approach to P. Now it seems to me that these two cases, which are quite 
distinct in theory, ought to manifest themselves in Nature and require distinguishing 
names. A race may be modified’ because a complex of organs with a certain system 
of values is good for it, or because it is bad for it. The race may be moditied because 
a certain element of it is fittest or because it is unfittest to survive. In the former 
case we select for survival round the centre, in the latter case we select for destruction. 
I propose to call these cases positive and negative selection respectively. It may be said 
that if there be positive selection in one part of the population there will be negative 
in another. But the kernel of the matter is in either case the existence of a centre, 
a definite set of most fit or of most unfit organs, while in positive selection the less fit 
organs, and in negative selection the more fit organs are distributed over wide areas 
of the field, and do not reach a maximum of unfitness or a maximum of fitness 
respectively for any definite individual. 

a 
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In the diagram we have also drawn the contour line to which the French population 
must extend if we are to get at least a representative population of 1000 Aino from 
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it, and further the contour line to which the Aino population must extend if we are 


to get at least a representative population of 1000 French from it. A consideration of 
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the nature of the contour lines of the surface of survivals shows that the contour lines 
above referred to, and marked “boundary” in the diagram, must touch the Aino 1 in 1000 
limit and the French 1 in 1000 limit respectively at the points in which they are 
touched by the contour lines p=1 and p’ = 1 of the corresponding surfaces of 
survivals. I have already indicated that the major axes of these boundaries are 
15:285 for the French and 15'890 for the Aino. The corresponding values of the 
parameter « are respectively given by 


15-285" | 15890 


Hence by (exviii.) we can easily find the frequency of population outside the 
contours kp and x,; 1f these be vp and v, we have : 


vy = (000,000,163, v, = '000,000,000,015. 


Thus the French population would have to be extended to a boundary in which only 
about 1 in six millions was excluded, and the Aino population to a boundary excluding 
only 15 in the billion! The boundaries of what we may thus term the selection 
populations are far larger than our conventional boundaries of 1 in 1000 for represen- 
tative populations. In fact, ıt would be impossible to select a representative Aino 
population from a conventional representative French population and vice versa—in 
either case the very exceptional members of French or Aino populations are 
required to complete the conventional representative populations of Aino or French 
by selection. 

(15.) I have devoted most of my consideration of the surface of survivals to 
a particular case in which two organs have been selected, and we consider the nature 
of p which determines the fraction of each group of individuals which survives. 
I have done this partly because normal surfaces are at best only an approximate 
representation of our selectable and selected distributions, and partly because I have 
thought a concrete case would best bring out the general points of investigations of 
this kind. 

But some little indication of the properties of the surface of survival-rates ought 
to be indicated here, or it may appear that they have been overlooked. While the 
contour lines of the correlation frequency surfaces for two organs must be ellipses, this 
does not follow in the case of the surface of survival-rates. In our illustration they 
were ellipses, but they may be also parabolas, hyperbolas, or even straight lines. 


H 


We must not therefore expect to find always a “centre” of positive or negative 


of contours with the rate 


H 


selection. We may come across a “saddle-back system ’ 
of survival constant along two intersecting lines, but rising in one pair of opposite 
angles and falling in the other pair. In this case we have fields of negative and 
positive selection separated by two independent relations between the two organs, 
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which are linear and for which the survival-rate is the same, they may be termed the 
“critical lines.” For one pair of angles the centre is now a “ centre of fitness,” for 
the other pair of angles a “centre of unfitness.” It seems to me that these critical 
organic relations may possess considerable biological importance. 

If the contour lines of the surface of survival-rates are parabolas, we have really 
only a limiting case of the centre at a very great distance. It is one in which the 
fittest (or most unfit) has no practical existence, but there is a direction towards 
which the rate of survival will be found to be always increasing or decreasing. 

If the contour lines of the surface of survival-rates are parallel straight lines, then 
so long as the deviation in one organ has a certain definite relation to that in the 
other, the survival-rate will remain constant. In this case the survival-rate will 
fall uniformly in one direction and remain constant in the direction at right angles 
to it. 

All the cases I have given here can occur just as easily as the elliptic contour 
system of our illustration and diagram. Each is marked by quite definite biological 
characteristics, and we may, perhaps, class them as elliptic, hyperbolic, parabolic, 
and linear selection. Even if the surface of survival-rates be not of the exponential 
quadric type discussed in this paper, yet to the neighbourhood of each part of it 
this classification of selection types will apply. 

If we pass to more than two organs, then similar considerations will apply; we 
shall only be reproducing the geometry of quadric surfaces in space of three and 
higher dimensions. But before we allow ourselves excursions into the higher 
geometry of the surface of survival-rates, it seems desirable that we should obtain 
quantitative determinations of this surface by experiments in artificial selection. 
We shall then be better able to see what part of our geometry will really be of 
service for the problems of natural selection. The field is too large to be cultivated 
for merely theoretical interests. We must first determine what parts of it are likely 
to have practical application to life as we find it, but of death-rates in the case of 
any living form but man, we are at present sadly ignorant. 
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